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SUMMARY OF RESULTS 


This program supported the analysis of IUE observations ot supernovae. One aspect was a 
Target-of -Opportunity program to observe bright supernovae which was applied to SN 1993J in M8 1 , 
and another was continuing analysis of the IUE data from SN 1987A. Because ot its quick response 
time, the IUE satellite has continued to provide useful data on the ultraviolet spectra of supernovae. 
Even after the launch of the Hubble Space Telescope, which has much more powerful ultraviolet 
spectrometers, the IUE has enabled us to obtain early and frequent measurements of ultraviolet 
radiation: this information has been folded in with our HST data to create unique observations of 
supernova which can be interpreted to give powerful constraints on the physical properties of the 
exploding stars. Our chief result in the present grant period was the completion of a detailed re- 
analysis of the data on the circumstellar shell of SN 1987 A. 

The presence of narrow high-temperature mission lines from nitrogen-rich gas close to SN 
1987 A has been the principal observational constraint on the evolution of the supernova s progenitor. 
Our new analysis shows that the onset of these lines, their rise to maximum, and their subsequent 
fading can be understood in the context of a model for the photoionization of circumstellar matter. 
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ABSTRACT 

Ultraviolet (X: 1 150-3300 A) observations of SN 1987A with the IUE satellite provide a unique data set. 
Observations started the day after discovery, 1 987 February 24 ( day 1.6), and a total of 75 1 spectra were obtained 
through 1992 June 9 (day 1567). The data have been processed to generate a complete catalog of SN 1987A 
ultraviolet spectra. The data reduction procedure includes careful line-by-line extraction, removing hits and hot 
pixels, and, most importantly, a scrupulous subtraction of the contribution from stars near SN 1 987 A within the 
IUE aperture. In addition to processing the data, we have also extracted light curves, and combined them with 
the ground-based optical data and II ST observations ( day 1278 to 243 1 ) to extend the study in both wavelength 
and time. The data-processing procedures of our IUE study produce results that are consistent with the HST data 
where they overlap, but not with the IUE study by Sanz Fernandez de Cordoba ( 1993) because of its incorrect 
background subtraction. The IUE data are consistent with the ground-based Walraven \ BLUW photometry, 
while we found that flux scale of the Soviet A STR O N satellite spectroscopic data is low by 15%. The UV flux 
plummeted during the earliest days of observations because of the drop in the photospheric temperature and the 
increase in opacity. However, after reaching a minimum of 0.04% on day 44, the UV flux increased by 1 75 times 
in its contribution to 7% of the total UVOIR bolometric luminosity at day 800. A revised set ot bolometric data 
has been constructed which includes the contribution of UV from day 1 to day 1352. Studies of the UV colors 
show that the supernova started to get bluer in UV around the time when dust started to form in the ejecta. Our 
results are consistent with the possibility that the dust condensed may be metal-rich. 

Subject headings: ultraviolet: stars — stars: individual (SN 1987A) 


t. INTRODUCTION 

SN I987A in the Large Magellanic Cloud (LMC) provided 
a unique opportunity for the study of supernovae which has 
resulted in unprecedented coverage in both wavelength and 
time. No other supernova has ever been as well observed as SN 
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1987A — it was truly a once-in-a-lifetime opportunity which 
we hope to repeat! SN 1987 A is the first supernova bright 
enough for long-term ultraviolet observations. Previous bright 
supernovae observed in the ultraviolet by the International Ul- 
traviolet Explorer (IUE), Type I Is SN 1979C and SN 1980K, 
and Type lb SN 1983N, only stayed bright enough for good 
signal-to-noise observations for less than a month ( Benvenuti 
et al. 1982). A large volume of literature has been generated by 
this spectacular event. Excellent reviews by Arnett et al. 
(1989), Imshennik & Nadezhin (1989), Hillebrandt & 
Hoflich ( 1989), McCray ( 1991 ), Nomoto et al. ( 1994), Che- 
valier (1992), Podsiadlowski (1992), and McCray ( 1993) 
cover various issues of the supernova, such as the neutrino 
burst, progenitor evolution, explosive nucleosynthesis, bolo- 
metric light curve, and interactions with circumstellar mate- 
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rial. In this paper we present the ultraviolet spectra of SN 
1987A taken by the IUE satellite from day 1 to day 1248, cov- 
ering the wavelength range 1 1 50-3300 A. The procedures used 
to process this unique data set are described. Beyond data pro- 
cessing, this set of data forms the basis for extensive analysis of 
the UV evolution of the supernova. 

Early results on the IUE spectroscopy were reported by 
Kirshner et al. ( 1987b) and Wamstaker et al. ( 1987b). These 
results stimulated a series of papers focussing on studying and 
interpreting the very early UV spectra of the supernova 
(Cassatella et al. 1987; Fransson et al. 1987; Panagia et al. 
1987b; Lucy 1987). The IUE spectra were also essential in 
identifying the supernova progenitor (Sonnebom, Altner, & 
Kirshner 1987; Gilmozzi et al. 1987). High-resolution IUE 
spectra (\/ AX » 10 4 ) with good signal-to-noise were obtained 
in the first few days of observations (Dupree et al. 1987; de 
Boer et al. 1987; Blades et al. 1988a,b; Savage et al. 1989). 
Interstellar absorption lines of various species with a wide 
range of ionization covering the heliocentric velocity range 
from -30 to 330 km s' 1 were detected. Absorptions originated 
from the galaxy and its halo {v < 120 km s _l ), LMC (u > 190 
km s' 1 ), and other sources other than the galaxy or LMC have 
been identified (Savage et al. 1989). Other studies of SN 
1987 A with the IUE satellite include the observations of the 
UV light echo in the vicinity of SN 1987A (Gilmozzi 1988; 
Panagia & Gilmozzi 1991 ). The UV echo spectrum could pro- 
vide information about the expansion velocity and the temper- 
ature of the supernova at the time of the shock reakout. 

Narrow emission lines of highly ionized species in the cir- 
cumsteliar matter were first detected in the UV spectra taken 
by the short-wavelength camera of IUE in 1987 May ( — day 
80)(Wamstekeretal. 1987a; Kirshner etal. 1987a; Sonnebom 
et al. 1995). High-resolution observation by IUE at day 275 
showed that the line widths of the C in] 1909, and the N III] 
1750 multiplets were ^30 km s' 1 ( Panagia et al. 1987a). Stud- 
ies of the relative strengths of the lines showed that the circum- 
stellar matter was highly CNO enriched ( Fransson et al. 1 989 ) . 
This is consistent with models where the progenitor was a red 
supergiant before exploding as a blue supergiant, thus losing 
much of its hydrogen envelope before the explosion. The emis- 
sion-line flux reached a maximum around day 410, and slowly 
decayed afterwards (Sonnebom 1991; Sonnebom et al. 1995; 
Sanz Fernandez de Cordoba 1993, hereafter SFdC93). Pa- 
nagia & Gilmozzi ( 1991 ) constructed a simple model for the 
emission-line light curves to determine the absolute diameter 
of the circular circumstellar ring. Combined with the measure- 
ment of the angular diameter by the Hubble Space Telescope 
(HST) (Jakobsen et al. 1991), the distance to SN 1987A, 
<4 n87a, can be determined (Panagia et al. 1991 ). The result, 
51.2 ± 3.1 kpc, is consistent with the other measurements of 
the distance to the LMC. Gould ( 1994a,b) adopted the emis- 
sion light curve model of Dwek & Felten ( 1 992 ) and obtained 
<^sn87a - 46.77 ± 0.76 kpc for a circular circumstellar ring and 
found that the results will be similar for an elliptical ring. De- 
tailed models of the circumstellar matter also exist ( Lundqvist 
1991; Dwek & Felten 1992; Plait et al. 1995) and they provide 
good agreement with the observed light curves of the IUE 
emission lines. The effects of these models on the determina- 
tion of ^sn8?a still remain to be seen. 

Summaries of early results from the IUE study of SN 1 987 A 
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can be found in Kirshner & Gilmozzi (1989), Gilmozzi 
( 1990), and Panagia & Gilmozzi ( 1991 ). A study similar to 
this one on the IUE data of SN 1987A can be found in the 
SFdC93 paper. In SFdC93, the IUE light curves of SN 1987A 
are computed from the time of the first IUE observation until 
1990 March (day 1100). The results and conclusions of the 
SFdC93 study are generally similar to this work. However, the 
UV flux attributed to SN 1987 A at late times in SFdC93 is 
about 7-15 times higher than what we find. The discrepancies 
are most probably caused by the improper subtraction of back- 
ground stars in SFdC93. The detailed differences between our 
data processing procedures and the ones employed in SFdC93 
are discussed in § 2. 

For this paper, we processed all the low-dispersion UV spec- 
tra of SN 1987A from the IUE archive. The data processing 
procedures we applied to all the 75 1 IUE spectra of SN 1 987 A 
are described in § 2. The resulting data set is significantly 
different from that in the IUE archive. The data have been 
carefully selected and combined to generate spectra with the 
best signal-to-noise. Moreover, as the supernova dims, remov- 
ing the flux from nearby field stars in the aperture becomes 
increasingly important. The extra flux comes mostly from Star 
2(V= 15.03) and Star 3 (V ~~ 15.7) in the Walker & Suntzeff 
( 1 990 ) catalog of stars near SN 1 987 A, which are always in the 
IUE aperture. Flux from Star 4, which is periodically present 
in the aperture as the spacecraft rotates annually around the 
sun, is also removed. It is only recently that we have obtained 
the definitive background spectra which permit this important 
correction to the data. The procedures used for removing these 
field stars are addressed in § 2. 

A total of 290 spectra spanning the epoch from day 1 to 1 248 
are assembled and are presented in § 3. This processed set of 
data forms the basis for our analysis of the U V evolution of the 
supernova. The UV spectra are combined with optical data to 
extend the wavelength coverage. We compare our data set to 
other UV observations of SN 1 987A from the Soviet ASTRON 
satellite (Lyubimkov 1990) and the Walraven VBLUW pho- 
tometry (Pel et al. 1987) in § 4. We show that the IUE results 
are consistent with the Walraven data within 5% and are higher 
than the ASTRON results by ~~ 1 5%. The UV light curves are 
shown in § 5 and are compared with the optical photometry. 
The ratio of ultraviolet flux to optical flux, F uv /F optical , is 
found to increase with time after about day 100. The pre- 
viously published UVOIR bolometric luminosity (Suntzeff & 
Bouchet 1990; Suntzeff et al. 1991; Bouchet, Danziger, & Lucy 
1991 ) includes the UV contribution for only the first week of 
IUE observations. We present the revised set of UVOIR bolo- 
metric luminosity, which includes the UV flux contributions 
from the IUE data from day 1 to 1352, in § 5. The UV color 
evolution is discussed in § 6. We found that the UV color be- 
came bluer just at the time that the dust started to condense in 
the SN 1987A ejecta. We conclude that the condensed dust 
does not resemble galactic or LMC interstellar dust in its ex- 
tinction properties and that it probably is metal-rich. Finally, 
in § 7, IUE results are compared with the HST Faint Object 
Camera observations ( day 1278-2431 ) obtained by the Super- 
nova INtensive Study (SINS) collaboration. We demonstrate 
that the IUE data archive can be connected to future UV ob- 
servations with imaging detectors for SN 1987 A and other 
supernovae. With the HST data, we establish a long-term rec- 
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ord of ultraviolet flux that will continue long after the super- 
nova became too faint for WE detections. Comparing the 
HST results with the WE data requires special attentions to 
the ‘Ted leaks” of the HST filters, which demand significant 
corrections to obtain valid results. We present a concluding 
discussion in § 8. 

2. UV OBSERVATIONS AND DATA REDUCTION 
2.1. WE Observations 

Taking advantage of a Target-of-Opportunity Program to 
observe bright supemovae, WE observations of SN 1987 A 
started at both NASA's Goddard Space Flight Center (GSFC) 
and ESA’s Villafranca del Castillo (VILSPA) promptly after 
its discovery. Subsequent systematic programs to obtain UV 
spectroscopy of SN 1987 A with WE have led to a thorough 
coverage of the ultraviolet evolution of the supernova. 

Detailed characteristics of the WE satellite, its scientific in- 
struments and performance can be found in Boggess et al. 

( 1978). One of the many attractive features of the geosynchro- 
nous WE satellite is that quick changes in observing schedules 
can be arranged by the ground observing stations which make 
the target-of-opportunity program very effective. The spectra 
are recorded on two image-converter SEC vidicon cameras: 
short wave prime (SWP) for the wavelength band 1 150-1975 
A, and long wave prime ( LWP ) for the wavelength band 1975- 
3300 A. We have processed all the spectra ofSN 1987A taken 
in the low dispersion mode ( resolution — 6-7 A ) with the large 
entrance aperture of size 10" X 20". In this mode, the spatial 
distribution of light projected on the long axis of the aperture 
is retained in the two-dimensional data. All WE spectra ob- 
tained at GSFC and VILSPA were processed by the WE Spec- 
tral Image Processing System ( IUESIPS) (Turnrose & Thomp- 
son 1984), as updated by Munoz Piero ( 1985) for extended 
line-by-line (ELBL) enhancements to improve the spatial res- 
olution of the low dispersion spectra. Each resulting ELBL file 
is a two-dimensional flux array with about 800 wavelength 
points (for SWP; 900 for LWP) and 110 lines in the spatial 
dimension, approximately parallel with the longer (20") axis 
of the aperture. The separation between spectra of two point 
sources in the ELBL file corresponds to the true angular sepa- 
ration in the sky projected on the long axis. The angular sepa- 
ration between two neighboring lines corresponds to 1 "07 ± 
0"04. In a typical large aperture observation of a point source, 
the spectral information is contained in about 18 lines in the 
ELBL file. Various spurious features in the low-dispersion WE 
spectra such as cosmic ray hits, hot pixels, and the reseaux 
etched on the faceplates of the cameras were flagged by the 
IUESIPS, enabling us to remove them from each individual 
line of the ELBL file where the supernova was present. This 
treatment is essential, especially in our late-time observations 
because of the diminishing flux of the supernova and the tong 
exposure time. 

Subsequent corrections are applied to adjust for the decreas- 
ing sensitivity with time of each WE camera. Teays & Garhart 
( 1990) showed that the LWP camera degradation in low dis- 
persion observing mode between the dates of 1984.5 and 
1989.5 can be represented by a linear decline. This is also true 
for the SWP camera degradation in low dispersion mode be- 
tween the dates of 1979.5 and 1992.0 (Garhart 1992a). Both 


studies indicated that the sensitivity changes varied signifi- 
cantly with wavelength due to the granular nature of the cam- 
era faceplate. Best-fit straight line coefficients were tabulated 
for the performance corrections at various wavelengths, and 
we have applied these corrections for all the WE spectra on SN 
1987 A, which spanned the period from 1987.2 to 1991.4. The 
latest broadband sensitivity monitoring confirmed the linear 
relation until year 1992.3 (Garhart 1992b) with a mean rate of 
degradation of —0.7%/year for the SWP camera and — 1.5%/ 
year for the LWP camera since the time of absolute calibration 
(1984.5-1985.0). 

2.2. Background Stars Subtraction 

SN 1987A is located in a densely populated region of the 
LMC and has close neighbors. Sonnebom et al. ( 1987 ) decon- 
volved the spatially resolved WE ELBL spectral files to sepa- 
rate the contributions from star 2 and star 3 present in the large 
aperture SWP data between day 30 and day 70. The absence of 
any hot stars left in the aperture confirmed the identification of 
star 1 of the Sanduleak -69°202 system as the progenitor that 
exploded as SN 1987A. Gilmozzi et al. (1987) conducted a 
similar analysis using the 1 300- 1 500 A data from day 7 to day 
14 and deduced the same progenitor. Moreover, Gilmozzi et 
al. ( 1987) observed stars 2 and 3 directly using the small WE 
aperture (diameter 3"). With star 2 located 2!9 ± 0"1 (PA 
315°) and star 3 located 1 ?6 ± 0T2 (PA 1 15°) from the super- 
nova, these two stars are always present in the large WE aper- 
ture ( 10" X 20"). Scuderi et al. ( 1995) compared the small 
aperture WE spectrum of star 2 with stellar atmosphere 
models with solar abundances and the WE spectra of standard 
stars, and they classified star 2 as a B1 III star. Walbom et al. 

( 1993) obtained the optical photometry and spectroscopy of 
stars 2 and 3. Star 2 is classified as a B2 III star, while star 3 
shows variations in UBVRI photometry as large as —0.3 mag, 
and in infrared photometry by about — 1 mag (Walbom et al. 
1993). Combining these results with the discovery of broad Fla 
emission in star 3 (Wang et al. 1992), it can be concluded that 
star 3 is a Be object. Be stars are known to vary in the UV as 
well ( Barylak & Doazan 1986), and thus star 3 complicates 
our study of the UV photometry of SN 1987A, as discussed in 
§5.1. 

Removing flux of the field stars in the spectra, which comes 
mostly from stars 2 and 3, grows more difficult and more im- 
portant as SN 1987 A fades. Contamination was especially se- 
vere in the SWP data where more than half of the integrated 
flux was from background stars as early as —20 days after the 
explosion. The WE spectra we obtained after day 1300 (for 
both the SWP and LWP cameras) remain constant in time ex- 
cept for the circumstellar emission lines, which suggests that 
they are completely dominated by flux from the background 
field stars instead of the supernova. These spectra (six in SWP 
and 1 1 in LWP between day 1 30 1 and 1 567 ) have been aver- 
aged, and carefully subtracted from all the data obtained. The 
combined background spectrum is shown in Figure 1. A de- 
finitive reduction of the WE data for SN 1987 A has not been 
possible without this measure of the background flux. 

Even after subtracting the background spectrum of Figure 1 
from earlier epochs, we still detected extra flux, with a six- 
month periodicity, in the remaining SWP spectra. When 



226 


PUN ET AL. 


Vol. 99 



Wavelength (A) 


Fig . 1 The average background spectru m ( day 1301-1567) subtracted from all the / UK data. The spectrum is not corrected for i nterstellar extinction 


pointed at SN I987A, the position angle of the WE aperture 
rotates once a year as the satellite is oriented for maximum 
power output from its solar cells. This is a natural way to pro- 
duce a 6 month periodic signal from additional field star(s) 
at a distance 5' to 10' from SN 1987A. We determined the 
maximum field star separation by studying the ELBL file and 
the position angle of the aperture at the epoch of maximum 
flux to locate the additional star(s) present. We conclude that 
the additional flux is from star 4 ( 575 from the supernova: V = 
18.24) in the Walker & Suntzeff ( 1990) catalog of stars near 
SN I987A. The extra flux from star 4 is removed by removing 
the lines in the ELBL spectral file which contribute additional 
flux and are spatially resolved from SN 1 987 A. star 2, and star 
3. All ELBL files of the SWP spectra taken after ~ day 400 are 
processed this way (and some before, at epochs of maximum 
star 4 contribution ) , when the fl ux from SN 1 98 7 A contributes 
only ~30% of the flux observed through the aperture and scru- 
pulous sky subtraction is required. 

Additional field stars located less than 5" from the supernova 
are always present in the WE aperture. Walker & Suntzeff 
( 1991 ) identified stars “4” and “5” located 2T0 and 2!9 away 
from the supernova. Notice that stars “4” and “5” in Walker 
& Suntzeff ( 1991 ) are different from stars 4 and 5 of the Walker 
& Suntzeff ( 1990) catalog of stars near SN 1987A. Plait et al. 
(1995) resolved star “5” into two stars with a separation of 
~0?1 in FOC images from HST. Also with the HST images. 
Plait, Chevalier, & Kirshner (1992) identified a star along our 
line of sight to the supernova superposed on the circumstellar 
ring located 0774 ± 0 TO 3 away. The colors and magnitudes are 
consistent with an A-type dwarf in the LMC ( Plait et al . 1 995 ) . 
Examining the HST images of SN 1987 A (described in § 7) 
shows that each of these stars is more than 5 mag dimmer than 
either star 2 or star 3 in the UV. Therefore, they are not sig- 
nificant sources of contamination to our data. 

The data processing procedures used in the SFdC93 paper 
are different from our present work in three major aspects. 


First, in SFdC93, not all the data from both the GSFC and 
V1LSPA stations between day 1 and 1 100 are included. Sec- 
ond, about 1 5% of the spectra included in the SFdC93 study to 
compute the U V light curves contained saturated pixels, which 
we have avoided in our data set. Third, and most important, we 
employed a different stellar background subtraction procedure 
from that of SFdC93. In SFdC93, by assuming the spectral 
types of stars 2 and 3 as B0 V and B 1.5 V, respectively, the 
combined background spectrum is generated with standard 
WE stars weighted according to the flux ratio deduced by Gil- 
mozzi et al. ( 1987). This combined spectrum is scaled to the 
UV flux in the late-time SN 1987 A spectra in the range 1300- 
1500 A and is subtracted from all the WE data as the back- 
ground. Because they assume a hotter atmosphere for star 2 
than others who have classified this star ( B2 III in Walbom et 
al. 1 993; B 1 III in Scuderi et al. 1 995 ) , the SFdC93 subtraction 
routine underestimates the long wavelength flux from the 
background stars. The star 2 and 3 background flux deter- 
mined by SFdC93 is about 1 1 % lower than we find in the range 
1 600- 1 950 A and about 25% lower in the range 1 952-3200 A. 
The measurement of the background flux is essential to the 
correct determination of the flux level from the supernova, es- 
pecially during late epochs when the background is larger than 
the supernova flux ( l > 6 days for SWP; t > 800 days for LWP) . 
At day 1000, the UV flux attributed to SN 1987 A by SFdC93 
is larger than ours by a factor of ~ 1 5 in the SWP range (1250- 
1 950 A ), and a factor of ~ 7 in LWP ( 2400-3200 A ) . 

2.3. Reddening Correction 

Because SN 1987A is only ~20' from the 30 Doradus neb- 
ula in the LMC, the supernova spectra have been corrected 
for interstellar extinction as derived for that region. Walker & 
Suntzeff( 1990) have measured the reddening in the field ofSN 
1 987 A by studying the U - B and B - V colors of 23 early-type 
stars in the 30 Doradus region. The reddening for SN 1 987A 
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was determined to be E( B - V ) total = 0.17 ± 0.02. Walborn et 
ai. ( 1989) studied the existing spectroscopic data of the pro- 
genitor of the supernova, Sk -69°202, and deduced a total red- 
dening toward SN 1987A of £*(2? ~ I r ) tota | = 0.16, with a galac- 
tic foreground reddening of E( B - V )gaiact.c — 0.06 and a LMC 
contribution of E(B - K) LMC ss 0.10. Gochermann, Goud- 
frooij, & Schmidt-Kaler ( 1989) have observed 4 1 galactic fore- 
ground stars in the direction of SN 1987A and deduced a ga- 
lactic foreground reddening of E{B — V )gaiact»c = 0.08 ± 0.01, 
which is higher than the usually assumed values. We have 
adopted the reddening values of E(B - E r ) galactic = 0.06 and 
E(B - L)lmc = 0.09 for direct comparisons with the results 
obtained by the CTIO group (Hamuy et al. 1988). However, 
by modeling the UV stellar atmosphere of star 2, Scuderi 
et al. ( 1995) deduced a total reddening of E(B - F r ) tota [ = 
0.19 ± 0.01, with —10% galactic contribution, i.e., 

FAB - f )galact»c = 0.02. 

The choice of dereddening function is important in UV 
studies because the correction is large and is known to vary 
in shape from place to place. Fitzpatrick ( 1985) documented 
significant differences between the ultraviolet extinction curves 
of the Galaxy, the LMC generally, and the 30 Doradus region. 
The 30 Doradus extinction curve has a weaker 2 1 75 A “bump" 
and a steeper far-ultraviolet rise than either the galactic or av- 
erage LMC extinction function. We used the 30 Doradus ex- 
tinction function evaluated by Fitzpatrick (1986) and assumed 
the ratio of visual extinction over the color excess, R v = A v j 
E(B - K), to be 3.1, the standard value for the diffuse inter- 
stellar medium. Fitzpatrick & Walborn ( 1 990 ) verified that the 
average 30 Doradus extinction correction function was justi- 
fied for Sk -69°203, a B0. 7 la supergiant, located 2'2 north of 
SN 1 987 A, by studying its UV and optical colors and spectra. 
This improves our confidence that the average 30 Doradus ex- 
tinction law is appropriate for the UV data of SN 1 987 A. For 
the galactic reddening component, the extinction function and 
the value of R y * 3.2, formulated by Seaton ( 1979) are used. 

3. UV SPECTROSCOPY OF SN 1987A 

3.1. UV Spectra 

We have processed the GSFC observations taken before 
1991 June 9 (day 1568), and VILSPA data up to 1991 March 
3 (day 1470). Some statistics of the low dispersion observa- 
tions we obtained from both GSFC and VILSPA are presented 
in Table 1. Because of the large wavelength-dependent varia- 
tions of sensitivities of both the SWP and LWP cameras 
(Harris & Sonneborn 1987), and the steep variations of flux 
with wavelength in the spectra, multiple exposures with differ- 

TABLE 1 


The Combined Statistics of the IUE Low-Dispersion SN 1 987 A 
Observations Obtained at GSFC and VILSPA through 1 992 June 



Number 
of SWP 
Spectra 

Number 
of LWP 
Spectra 

Total 

IUE 

Spectra 

Observation total 

221 

530 

751 

Background (> day 1 300) 

6 

1 1 

17 

Raw data (< day 1300) 

215 

519 

734 

Combined data (day 1-1248) 

124 

166 

290 


ent exposure times were often needed. We have combined 
them to ensure optimum coverage of a wavelength band with 
satisfactory signal-to-noise ratio. However, in order to preserve 
information on the time evolution of the supernova, only ob- 
servations made within 5 hr were considered for combination. 
Thus many of the processed IUE spectra are combined from 
two or, less often, three observations. There were only two oc- 
casions when we combined more than three spectra. 

The catalog of the SWP spectra taken between day 1 30 1 and 
day 1568 which we used to create the background spectrum 
is given in Table 2A. The corresponding catalog for the LWP 
spectra is presented in Table 2B. The third column is the time 
of observation in days after core collapse of the supernova, 
1987 February 23.316 (UT)(JD = 2,446,849.816), as deter- 
mined at the IMB and Kamiokande II neutrino detectors 
(Bionta et al. 1987; Hirata et al. 1987). Complete lists of IUE 
observations of the SWP and LWP cameras on SN 198 7 A are 
given in Table 3A and 3B, respectively. The exposure time in 
the fourth column is corrected for the IUE exposure time 
quantization and the camera response time (Crenshaw 1988). 
These corrections are important for the early data where the 
time of exposure is short (t c%p < 1 5 s). A total of 290 combined 
spectra are obtained from day 1 to day 1248. The IUE spectra 
at different epochs are presented in Figure 2, where the number 
to the right of a spectrum represents the number of days after 
outburst and the displayed spectrum is the average over the 
indicated time period. The vertical axis is the logarithm of ob- 
served flux and successive spectra are shifted vertically for clar- 
ity of display, except for the spectra taken in the first four days. 
The magnitude of shift, in dex, is shown in the parenthesis un- 
der the epoch of observation. 

There had been five Type II supemovae (SNe II) observed 
by IUE before SN 1987 A. However, only two, SN 1979C and 
SN 1980K. were bright enough for high-quality UV data to 


TABLE 2 

Spectra Averaged for Background Subtraction 


IUE Image 

Date (UT) 

Day 

^cxp 

(minutes) 

Short-wavelength spectra ( 1 1 50- 1 975 A) 

SWP 39757 

1990 Oct 4.02 

1318.71 

265 

SWP 40002 

1990 Oct 31.03 

1345.72 

270 

SWP 40275 

1990 Dec 5.85 

1381.54 

280 

SWP 40858 

1991 Feb 12.63 

1450.32 

280 

SWP41 179 

1991 Mar 25.68 

1491.37 

310 

SWP 4 1802 

1991 Jun 9.34 

1567.03 

290 


Long-wavelength spectra ( 1 975-3300 A) 


LWP 18792 

1990 Sep 15.57 

1300.26 

110 

LWP 18933 

1990 Oct 4. 16 

1318.85 

120 

LWP 19090 

1990 Oct 31.17 

1345.86 

105 

LWP 19356 

1990 Dec 5.99 

1381.68 

120 

LWP 19518 

1991 Jan 7.42 

1414.11 

130 

LWP 19743 

1991 Feb 12.77 

1450.46 

120 

LWP 19865 

1991 Mar 3.21 

1468.90 

130 

LWP 19985 

1991 Mar 25.52 

1491.21 

140 

LWP 20 100 

1991 Apr 7.95 

1504.64 

90 

LWP 20 102 

1991 Apr 8.49 

1505.18 

140 

LWP 20553 

1991 Jun 9.48 

1567.17 

125 


TABLE 3 

WE Camera Observations of SN 1 987 A 


WE Image 


Date (UT) Day 

SWP camera (1150-1 975 A) 


SWP 30376 

SWP 30378 

SWP 30380 

SWP 30385 

SWP 30388 

SWP 30390 

SWP 30395 

SWP 30397 

SWP 30398 

SWP 30401,30402 
SWP 30407, 30408 
SWP 30410, 30411 
SWP 30414, 30415 
SWP 30416, 30417 

SWP 30421 

SWP 30422 

SWP 30427 

SWP 30428 

SWP 30429 

SWP 30433 

SWP 30472 

SWP 305 12 

SWP 30522 

SWP 30547 

SWP 30551 

SWP 30592 

SWP 30637 

SWP 30743 

SWP 30841 

SWP 30871 

SWP 30882 

SWP 30907 

SWP 30929 

SWP 30932 

SWP 30974 

SWP 3 1000 

SWP 3 1040 

SWP 3 1064 

SWP 31 125 

SWP 31132 

SWP 31 154 

SWP 3 1166 

SWP 31 177 

SWP 3 1273. 31274 . 
SWP31319, 31320 . 

SWP 3 1334 

SWP31371, 31372 . 
SWP 3 1420, 31421 . 
SWP 3 1462, 31463 . 
SWP 31818, 31819 . 
SWP 32030, 32031 . 
SWP 322 19, 32220 . 

SWP 32395 

SWP 32619, 32620 . 
SWP 32797, 32798 . 

SWP 32879 

SWP 329 10, 329 1 1 . 

SWP 32984 

SWP 33035, 33036 . 
SWP 33 104, 33105 . 
SWP 33 175, 33 176 . 
SWP 33279, 33280 . 
SWP 33331, 33332 . 

SWP 33423 

SWP 33492 

SWP 33496, 33497 .. 
SWP 335 19, 33520 . 


1987 Feb 24.87 

1.55 

10 s 

1987 Feb 24.98 

1.67 

8 

1987 Feb 25. 16 

1.84 

10 

1987 Feb 25.44 

2.12 

10 

1987 Feb 25.86 

2.54 

10 

1987 Feb 25.95 

2.63 

18 

1987 Feb 26.23 

2.92 

20 

1987 Feb 26.43 

3.11 

30 

1987 Feb 26.51 

3.19 

30 

1987 Feb 26.84 

3.53 

1, 25 minutes 

1987 Feb 27.53 

4.22 

4.5, 90 

1987 Feb 27.95 

4.63 

10,90 

1987 Mar 1.01 

5.69 

90,40 

1987 Mar 1.28 

5.97 

40, 170 

1987 Mar 2.02 

6.70 

90 

1987 Mar 2.25 

6.93 

106 

1987 Mar 3.31 

8.00 

161 

1987 Mar 3.88 

8.56 

180 

1987 Mar 4.00 

8.68 

198 

1987 Mar 5.38 

10.06 

195 

1987 Mar 9.31 

13.99 

180 

1987 Mar 13.59 

18.27 

240 

1987 Mar 14.37 

19.05 

182 

1987 Mar 16.61 

21.29 

240 

1987 Mar 17.24 

21.93 

87 

1987 Mar 22.53 

27.22 

300 

1987 Mar 28.37 

33.05 

185 

1987 Apr 8.55 

44.23 

225 

1987 Apr 22.76 

58.44 

45 

1987 Apr 26.62 

62.30 

40 

1987 Apr 29.79 

65.48 

110 

1987 May 4.56 

70.25 

230 

1987 May 7.70 

73.39 

60 

1987 May 8.10 

73.78 

60 

1987 May 14.91 

80.59 

120 

1987 May 19.54 

85.23 

180 

1987 May 25.17 

90.86 

195 

1987 May 30.57 

96.26 

180 

1987 Jun 9.08 

105.77 

200 

1987 Jun 11.13 

107.81 

195 

1987 Jun 14. 14 

110.82 

180 

1987 Jun 16.30 

1 12.99 

180 

1987 Jun 17.07 

113.76 

128 

1987 Jul 1.38 

128.06 

240, 140 

1987 Jul 1 1.99 

138.68 

240, 120 

1987 Jul 12.74 

1 39.42 

80 

1987 Jul 20.35 

147.03 

240, 105 

1987 Jul 29.06 

154.75 

240, 120 

1987 Aug 3.27 

160.95 

240, 85 

1987 Sep 10.12 

198.81 

90, 240 

1987 Oct 9.07 

227.75 

240, 90 

1987 Nov 3.02 

252.70 

240, 80 

1987 Nov 26.08 

275.76 

90 

1987 Dec 25.92 

305.61 

240, 80 

1988 Jan 27.83 

338.51 

210, 70 

1988 Feb 7.81 

349.49 

75 

1988 Feb 13.32 

355.00 

70,212 

1988 Feb 25.66 

367.34 

60 

1988 Mar 4.68 

375.37 

240, 40 

1988 Mar 17.27 

387.96 

70, 240 

1988 Mar 29.66 

400.34 

210, 90 

1988 Apr 14.24 

415.92 

60, 200 

1988 Apr 2 1.59 

423.28 

240,75 

1988 May 1.91 

433.59 

80 

1988 May 8.92 

440.60 

80 

1988 May 10.11 

441.79 

70, 240 

1988 May 13.49 

445.18 

240, 80 


228 


TABLE 3 — Continued 


WE Image Date(UT) Day 

SWP 33644 1988 May 26.61 458.29 80 

SWP 33725 1988 Jun 8.74 471.42 80 

SWP 3374 1, 33742 1988 Jun 11.02 473.71 80,235 

SWP 33799, 33800 1988 Jun 22.40 485.09 240,80 

SWP 33832 1988 Jun 29.65 492.33 65 

SWP 33868, 33869 1988 Jut 5.95 498.64 90, 195 

SWP 33899 1988 Jul 12.55 505.24 80 

SWP 33937. 33938 1988 Jul 18.32 511.00 240,80 

SWP 33966, 33967 1988 Jul 25.95 518.63 90,251 

SWP 34057, 34058 l988 Aug9.24 532.93 240,80 

SWP 34087, 34088 1988 Augl4.73 538.42 60,200 

SWP 34095 1988 Aug 17.34 541.02 110 

SWP 34193 1988 Sep 7.26 561.95 120 

SWP 34231, 34232 1988 Sep 12.77 567.46 80,200 

SWP 34316, 34317* 1988 Sep 26. 19 580.87 120,110 

SWP 34441 1988 Oct 8.81 593.50 150 

SWP34477 I9880ctl5.47 600.16 110 

SWP 34670 1988 Nov 3.58 619.26 133 

SWP 34701 1988 Nov 7.40 623.09 120 

SWP 34751, 34752 1988 Novl7.05 632.73 240,42 

SWP 34807 1988 Nov 24.40 640.09 135 

SWP 3487 1 1988 Dec 4.98 650.66 163 

SWP 35030 1988 Dec 16.00 660.68 110 

SWP35096 1988Dec22.32 668.01 130 

SWP 35 126 1988 Dec 24.97 670.66 180 

SWP 35243 1989 Jan 4.05 680.74 90 

SWP 35308 1989 Jan 10.87 687.55 220 

SWP 35379 1989 Jan 19.06 695.75 140 

SWP 35401 1989 Jan 24.53 701.21 160 

SWP 35505 1989 Feb 7.59 715.28 220 

SWP35686 1989 Mar5.61 741.29 240 

SWP 35822 1989 Mar 20.39 756.07 183 

SWP 35940 1989 Apr 5.50 772.18 260 

SWP36035 1 989 Apr 19.3 1 785.99 166 

SWP 36171 1989 Apr 30.47 797.16 210 

SWP 36258 1989 May 14.21 810.89 222 

SWP36279 l989May 17.47 814.15 260 

SWP 36539 1989 Jun 17.34 845.03 260 

SWP 36578 1989 Jun 26.09 854.78 310 

SWP36670 1989 Jul 13.39 871.07 200 

SWP 36676 1989 Jul 14.26 871.95 260 

SWP 36801 1 989 Aug 6. 18 894.87 260 

SWP 36968 1989 Sep 9.85 929.54 270 

SWP 37062 1989 Sep 17.80 937.49 230 

SWP 37088 1989 Sep 20.12 939.80 300 

SWP 37236 1989 Oct 2.78 952.46 263 

SWP 37424 1989 Oct 22.00 971.69 230 

SWP 37574 1989 Nov 12.95 993.63 300 

SWP37798 1989 Dec 11.87 1022.55 300 

SWP 37973 1990 Jan 8.77 1050.46 285 

SWP38055 1990Jan20.45 1062.13 280 

SWP 38172 1990 Feb9.62 1082.30 285 

SWP 38307 1990 Mar 5.61 1106.29 285 

SWP 38336 1990 Mar 12.35 1113.04 275 

SWP 38536 1990 Apr 6.52 1138.20 285 

SWP 38555 1 990 Apr 9.27 1140.96 270 

SWP 38866 1990 May 25.44 1187.13 285 

SWP 39300 1990 Jul 25.29 1247.97 285 

LW P camera ( 1 975-3300 A) 

LWP I0189 a 1987 Feb 24.88 1.56 1.9s 

LWP 10191 1987 Feb 24.97 1.66 Lis 

LWP 10193, 10195, 10198 1987 Feb 25.28 1.96 1.1s, 3.6s, 1- Is 

LWP 10199 1987 Feb 25.83 2.51 Lis 

LWP 10202, 10203, 10205 l987 Feb26.26 2.95 Lis, 1 3s, 2.7s 

LWP 10207 1987Feb26.50 3.18 2.3s 

LWP 10210, 1021 1 1987 Feb 26.83 3.52 2.7s, 90s 

LWP 10218 1987 Feb 27.50 4.18 5.6s 

LWP 10220, 10221, 10222 1987 Feb 27.84 4.53 3m, 7.7s, 4m 
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TABLE 3 — Continued 


WE Image 


LWP 10227, 10228 

LWP 10230, 10231 

LWP 10233, 10234 

LWP 10241, 10242, 10243, 10244 

LWP 10248, 10249 

LWP 10251, 10252, 10253, 10254, 10255 

LWP 10259, 10260 

LWP 10270, 10271, 10272 

LWP 10287, 10288, 10289 

LWP 10300, 10301 

LWP 10302, 10303, 10304 

LWP 10318, 10319 

LWP 10322, 10323, 10324 

LWP 10344, 10345, 10346 

LWP 10349, 10350 

LWP 10351, 10352 

LWP 10371, 10372 

LWP 10388, 10389, 10390 

LWP 10401, 10402 

LWP 10441, 10442 

LWP 10445, 10446 

LWP 10518, 10519, 10520 

LWP 10587, 10588 

LWP 10647, 10648 

LWP 10664, 10665 

LWP 10669, 10670 

LWP 10683, 10684 

LWP 10692, 10693 

LWP 10706, 10707 

LWP 10711, 10712 

LWP 10723, 10724 

LWP 10755, 10756 

LWP 10765, 10766 

LWP 10788, 10789 

LWP 10811, 10812 

LWP 10826, 10828 

LWP 10833, 10834 

LWP 10873, 10874, 10875 .. . 

LWP 10888, 10889 

LWP 10909, 10910 

LWP 10939, 10940 

LWP 10960, 10961 

LWP 10970, 10971 

LWP 11080, 11081 

LWP 11086, 11087 

LWP 11130, 11131 

LWP 11 182, 11183 

LWP 11201, 11202 

LWP 11225. 11226 

LWP 11254, 11255 

LWP 11281, 1 1282 

LWP 11311, 11312 

LWP 11 369, 11370 

LWP 11436, 11437 

LWP 1 1500, 1 1501, 1 1502 ... 

LWP 11503, 11504 

LWP 1 1534, 11535 

LWP 11588, 11589 

LWP 11672, 11673 

LWP 11824, 11825, 11826 

LWP 11953, 11954 

LWP 11997, 11998 

LWP 12083, 12084 

LWP 12167, 12168 

LWP 12235, 12236 

LWP 12294, 12295 

LWP 12372, 12373 

LWP 12410. 12411 

LWP 12486, 12487 

LWP 12525, 12526 


Date(UT) Day 


1987 Feb 28.96 

5.65 

1987 Mar 1.27 

5.95 

1987 Mar 1.37 

6.06 

1987 Mar 2.24 

6.92 

1987 Mar 2.84 

7.53 

1987 Mar 3.25 

7.93 

1987 Mar 4.01 

8.69 

1987 Mar 4.91 

9.60 

1987 Mar 6.72 

11.41 

1987 Mar 8.74 

13.43 

1987 Mar 9.33 

14.02 

1987 Mar 13.53 

18.21 

1987 Mar 14.26 

18.94 

1987 Mar 16.58 

21.26 

1987 Mar 17.18 

21.86 

1987 Mar 17.30 

21.99 

1987 Mar 19.84 

24.52 

1987 Mar 22.59 

27.27 

1987 Mar 23.85 

28.53 

1987 Mar 28.33 

33.02 

1987 Mar 28.97 

33.66 

1987 Apr 8.54 

44.22 

1987 Apr 16.66 

52.35 

1987 Apr 26.66 

62.33 

1987 Apr 29.80 

65.48 

1987 Apr 30.73 

66.42 

1987 May 2.84 

68.53 

1987 May 4.61 

70.30 

1987 May 7.71 

73.39 

1987 May 8.03 

73.71 

1987 May 9.67 

75.35 

1987 May 13.63 

79.32 

1987 May 14.86 

80.55 

1987 May 19.52 

85.20 

1987 May 22.68 

88.36 

1987 May 25.13 

90.81 

1987 May 25.82 

91.51 

1987 May 30.39 

96.08 

1987 Jun 1.73 

98.41 

1987 Jun 4.74 

101.42 

1987 Jun 7.60 

104.29 

1987 Jun 9.10 

105.79 

1987 Jun 10.94 

107.63 

1987 Jun 25.05 

121.73 

1987 Jun 25.78 

122.46 

1987 Jul 1.34 

128.03 

1987 Jul 10.94 

137.62 

1987 Jul 12.74 

139.42 

1987 Jul 17.01 

143.70 

1987 Jul 20.27 

146.96 

1987 Jul 26.95 

153.63 

1987 Aug 3. 14 

160.83 

1987 Aug 10.84 

168.53 

1987 Aug 20.04 

177.73 

1987 Aug 27.74 

185.43 

1987 Aug 27.96 

185.64 

1987 Sep 1.17 

188.86 

1987 Sep 10.05 

198.74 

1987 Sep 20.74 

209.42 

1987 Oct 9.06 

227.75 

1987 Oct 26. 15 

244.83 

1987 Nov 2.89 

252.58 

1987 Nov 13.68 

263.37 

1987 Nov 25.98 

275.66 

1987 Dec 6. 18 

285.87 

1987 Dec 16.49 

296.18 

1987 Dec 25.86 

305.54 

1988 Jan 2.03 

312.71 

1988 Jan 13.54 

324.22 

1988 Jan 19.13 

329.82 


tex p 

28s, 12m 

lm 16s, 4m 

6m 30s, 1 m 1 6s 

2m 20s, 2m 30s, 2m 30s, 4m 

2 m 50s, 10m 

lm 25s, 2m, 90m, 2m 10s, 2m 30s 

40m, 5m 

6m, 15m, 60m 

4m 20s, 120m, 18m 

6m, 50m 

7m, 6m, 100m 

7m, 24m 

1 5m, 7m, 90m 

7m, 24m, 1 10m 

6m 45s, 50m 

6m 55s, 90m 

7m, 24m 

7m, 24m, 90m 

7m, 24m 

230m, 6m 

6m, 20m 

5m 30s, 20m, 1 1 5m 
5m 30s, 120m 
4m, 55m 
3m 35s, 55m 
3m 35s, 25m 
3m, 25m 
2m 45s, 22m 
2m 30s, 20m 
3m 30s, 90m 
2m 30s, 20m 
2m 20s, 17m 
2m 20s, 20m 
2 m 20s, 20m 
2 m 20s, 20m 
2m 30s, 60m 
2m 20s, 15m 
6m 30s, 39m, 120m 
2m 35s, 15m 
2m 40s, 13m 
2m 40s, 1 3m 
2m 30s, 60m 
2m 50s, 17m 
4m 10s, 1 10m 
4m 5s, 20m 
4m 10s, 20m 
4m 20s, 25m 
4m 20s, 20m 
20m, 4m 
3m 50s, 16m 
4m 30s, 35m 
3m 40s, 15m 
4m 20s, 40m 
3m 30s, 13m 
3m 40s, 20m, 90m 
3m 10s, 12m 
3m 10s, 12m 
3m 10s, 12m 
3m, 20m 
3m, 12m, 35m 
3m, 12m 
3m, 12m 
2m 45s, 20m 
2m 50s, 10m 
3m, 12m 
2m 30s, 10m 
3m, 12m 
3m, 12m 
3m 30s, 1 2m 
12m, 3m 
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TABLE 3 — Continued 


WE Image 

Date (UT) 

Day 


LWP 12558, 12559 

1988 Jan 24.17 

334.85 

3m, 12m 

LWP 12576, 12577 

1988 Jan 27.79 

338.47 

3m, 60m 

LWP 12623, 12624 

1988 Feb 7.79 

349.48 

3m, 12m 

LWP 12656, 12657 

1988 Feb 13.23 

354.92 

3m, 25m 

LWP 12677, 12678 

1988 Feb 17.03 

358.72 

3m, 12m 

LWP 12734, 12735 

1988 Feb 25.73 

367.41 

1 2m. 3m 

LWP 12789, 12790 

1988 Mar 3.21 

373.89 

3m, 1 2m 

LWP 12869, 12870 

1988 Mar 17.22 

387.90 

3m, 30m 

LWP 12942, 12943, 12944 

1988 Mar 29.68 

400.36 

3m. 12m, 30m 

LWP 12991, 12992 

1988 Apr 5,96 

407.65 

3m 30s, 12m 

LWP 13041, 13042 

1988 Apr 14,19 

415.87 

3m 30s, 22m 

LWP 13087, 13088 

1988 Apr 21.57 

423.26 

3m 30s, 12m 

LWP 13153, 13154 

1988 May 1.92 

433.61 

4m, 12m 

LWP 13193, 13194, 13195 

1988 May 10.06 

441.74 

4m, 30m, 20m 

LWP 13217, 13218 

1988 May 13.42 

445.11 

4m. 12m 

LWP 13236, 13237 

1988 May 15.88 

447.57 

4m, 1 2m 

LWP 13310, 13311 

1988 May 26.58 

458.27 

12m, 4m 

LWP 13387, 13388 

1988 Jun 8.74 

471.42 

4m, 12m 

LWP 13402, 13403, 13404 

1988 Jun 11.10 

473.79 

4m, 60m, 20m 

LWP 13489, 13490 

1988 Jun 22.34 

485.02 

4m, 12m 

LWP 13528, 13529 

1988 Jun 29,62 

492.31 

60m, 4m 

LWP 13573, 13574, 13575 

1988 Jul 5.96 

498.64 

4m 30s, 60m, 23m 

LWP 13624. 13625 

1988 Jul 12.55 

505.24 

4m, 15m 

LWP 13684, 13685 

1988 Jul 1 8.25 

510.94 

4m 30s, 1 3m 

LWP 13741. 13742 

1988 Jul 25.88 

518.56 

4m 30s, 40m 

LWP 13743, 13744 

1988 Jul 26.52 

519.20 

4m 30s, 13m 

LWP 13787, 13788 

1988 Aug 4.66 

528.34 

5m, 15m 

LWP 13813, 13814 

1988 Aug 8. 18 

532.86 

5m, 15m 

LWP 13846, 13847, 13848 

1988 Aug 14.75 

538.43 

4m 30s, 30m, 55m 

LWP 13866, 13867 

1988 Aug 17.34 

541.03 

6m, 18m 

LWP 13896, 13897, 13898 

1988 Aug 22.54 

546.22 

6m, 18m, 60m 

LWP 14007, 14008, 14009 

1988 Sep 7.35 

562.04 

6m 30s, 18m, 80m 

LWP 14038, 14039, 14040 

1988 Sep 12.79 

567.47 

4m 30s, 25m, 45m 

LWP 14126, 14127 

1988 Sep 26. 14 

580.82 

6m 30s, 18m 

LWP 14190, 14191 

1988 Oct 8.71 

593.39 

8m, 25m 

LWP 14229, 14230 

1988 Oct 15.45 

600.13 

7m 30s, 25m 

LWP 14348, 14349 

1988 Oct 31.08 

615.76 

12m, 60m 

LWP 14476, 14477 

1988 Nov 17.04 

632.72 

9m, 25m 

LWP 14692 

1988 Dec 22.26 

667.94 

12m 

LWP 14766 

1989 Jan 4.06 

680.74 

I4m 

LWP 14799 

1989 Jan 10.69 

687.37 

1 5m 30s 

LWP 14865 

1989 Jan 19.00 

695.69 

16m 

LWP 14909, 14910 

1989 Jan 24.53 

701.22 

13m, 45m 

LWP 14981, 14982 

1989 Feb 7.71 

715.40 

19m, 80m 

LWP 15059, 15060 

1989 Feb 19.01 

726.70 

21m, 53m 

LWP 15133, 15134 

1989 Mar 5.73 

741.42 

27m, 65m 

LWP 15191, 15192 

1989 Mar 13.88 

749.56 

3 1 m, 90m 

LWP 15216, 15217 

1989 Mar 18.89 

754.58 

32m, 80m 

LWP 15225, 15226 

1989 Mar 20.25 

755.93 

60m, 15m 

LWP 15307, 15308 

1989 Apr 5.63 

772.32 

35m, 35m 

LWP 15375 

1989 Apr 19.23 

785.92 

40m 

LWP 15391, 15392 

1989 Apr 22.79 

789.48 

42m, 1 1 5m 

LWP 15414, 15415 

1989 Apr 30.60 

797.28 

40m, 90m 

LWP 15512 

1989 May 14.10 

810.79 

40m 

LWP 15530, 15531 

1989 May 17.48 

814.16 

45m, 78m 

LWP 15641 

1989 Jun 2.62 

830.30 

50m 

LWP 15741 

1989 Jun 17.46 

845.14 

55m 

LWP 15771, 15772 

1989 Jun 22.80 

850.48 

55m, 95m 

LWP 15796 

1989 Jun 26.96 

854.64 

60m 

LWP 15855, 15856 

1989 Jul 5.82 

863.41 

55m, 75m 

LWP 15914, 15915 

1989 Jul 14.41 

872.09 

55m, 46m 

LWP 16076 

1989 Aug 6.30 

894.99 

70m 

LWP 16185 

1989 Aug 2 1.38 

910.06 

150m 

LWP 16235 

1989 Aug 28.52 

917.21 

80m 

LWP 16295, 16296 

1989 Sep 7.55 

927.23 

80m, 80m 

LWP 16313 

1989 Sep 9.72 

929.41 

90m 

LWP 16372 

1989 Sep 17.69 

937.37 

70m 

LWP 16380 

1989 Sep 20.26 

939.94 

85m 

LWP 16466 

1989 Oct 2.65 

952.33 

90m 

LWP 16527 

1989 Oct 10.44 

960.12 

95m 
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TABLE 3 — Continued 


WE Image Date(UT) Day /„ p 


LWP 16600 1989 Oct 22.12 971.80 80m 

LWP16751 1989 Nov 6.22 986.90 110m 

LWP 16780 1989 Nov 13.09 993.77 95m 

LWP 16841 1989 Nov 26.39 1007.07 100m 

LWP 16855 1989 Nov 28.54 1009.22 100m 

LWP 16911 1989 Dec 12.01 1022.69 90m 

LWP 171 13 1990 Jan 8.91 1050.60 110m 

LWP 17189 1990 Jan 20.58 1062.27 94m 

LWP 17329 1990 Feb 9.76 1082.44 105m 

LWP 17413 1990Feb 23.65 1096.33 120m 

LWP 17475 1990 Mar 5.75 1106.43 105m 

LWP 17519 1990 Mar 12.22 1112.91 90m 

LWP 17543 1990 Mar 17.07 1117.75 150m 

LWP 17702 1990 Apr 6.66 1138.35 112m 

LWP 17714 1990 Apr 9. 13 1140.82 120m 

LWP 17783 1990 Apr 21.99 1153.67 140m 

LWP 17987 1990 May 25.58 1187.27 100m 

LWP 18440 1 990 Jul 25.42 1248.11 85m 


* Spectrum saturated for 2650-3000 A. 


be taken for line identification and spectral evolution studies 
(Benvenuti et al. 1982). During the first two weeks of IUE 
observations, the UV spectra from both these supemovae 
showed a strong UV continuum below 1 500 A with an excess 
over blackbody radiation extrapolated from the optical spec- 
tra. Furthermore, strong emission lines of highly ionized spe- 
cies such as N v 1240, Si iv 1400, and N in] 1750 were ob- 
served blueward of ~2000 A. Fransson et al. ( 1984) showed 
that the UV excess may originate from inverse Compton scat- 
tering by energetic, thermal electrons ( T ~~ 10 9 K) produced 
just behind the outward moving shock front where the super- 
nova ejecta interacts with the preexisting circumstellar enve- 
lope. The radiation from this shock wave dominates the con- 
tinuum below 1500 A, while the Comptonized flux which 
scatters back to the interior is the dominant source of energy 
for the UV emission lines. These high-energy photons cause 
the ionization and excitation of the species in the outer layers 
of the supernova envelope, where the density is lower and the 
ionizing radiation flux is higher than in the inner envelope, 
producing narrow emission lines through fluorescence. 

The UV spectra from SN 1987A showed marked contrasts 
with these two “classical” SNe II. A large UV flux was seen in 
the first IUE spectra obtained (day 1.6) with a high photo- 
spheric temperature of about 14,000 K. The broad P Cygni 
features of Mg n at ~~2800 A indicated high expansion veloci- 
ties in the debris extending up to 40,000 km s “ 1 . The extra U V 
continuum below 1 500 A and the highly ionized UV emission 
lines blueward of 2000 A are missing in the spectra of SN 
1 987 A, suggesting less interaction of the supernova ejecta with 
the intermediate circumstellar environment than SN 1979C 
and SN 1 980K. In the next few days, a huge drop in flux is seen 
at all UV wavelengths. In less than 3 days (day 4.6), radiation 
blueward of 1600 A decreased by a factor of 3000. At the same 
time, the spectral features changed dramatically: instead of a 
continuum, the spectrum consisted of broad absorption bands 
caused by iron peak elements. The spectral shape is very sim- 
ilar to that of Type la supemovae ( SNe la), such as SN 1980N, 
and SN 198 1 B ( Panagia & Gilmozzi 1991). Ultraviolet spectra 
were obtained for SN 1992A, another SN la, by both IUE 


(^2000-3300 A) and the Faint Object Spectrograph (FOS) 
on HST ( ~ 1800-4800 A) (Kirshner et al. 1993). The HST 
UV spectrum of SN 1 992A taken 5 days past optical maximum 
is strikingly similar to that of SN 1 987A observed at day 3.5 in 
the range 1800-3200 A (Jeffery et al. 1994, Fig. 3), with both 
spectra showing broad absorptions at around 1930, 2180, 
2380, 2600, and 3050 A. In contrast, the UV spectrum for SN 
1993J, a SN lib, taken by FOS of HST about 18 days after 
the explosion does not show any obvious broad absorptions 
(Jeffery et al. 1994), and instead resembles SN 1979C and SN 
1980K. Earlier UV spectra of SN 1993J taken with IUE also 
lack obvious broad absorptions and show only a strong UV 
continuum (Wamsteker et al. 1993; Sonnebom et al. 1993). 
These observations led Jeffery et al. ( 1994) to suggest that the 
interaction of the supernova ejecta with the circumstellar 
matter surrounding the progenitor is the source of differences 
in the UV behavior of supemovae. A smooth UV spectrum 
lacking in broadband absorptions caused by iron group ele- 
ments is characteristic of supemovae which interacted with 
thick circumstellar envelopes. In the case of SN 1987 A, the 
circumstellar material is present in the form of a ring, located 
200 light-days from the supernova (Jakobsen et al. 1991; Plait 
et al. 1995), with little material in the immediate neighbor- 
hood. The ejecta of SN 1987 A had relatively little interaction 
with the circumstellar matter compared to SN 1993J, SN 
1 980K, or SN 1 979C. This may explain why the U V spectrum 
of SN 1987A resembles SNe la more than other SNe II. 

The UV opacity is primarily that of line opacity. There are 
many strong and closely spaced resonance lines in the UV re- 
gion, especially from the neutral or singly ionized iron group 
elements, such as Fe, Ni, and Ti. Assuming LTE, Branch 
( 1987) computed the Sobolev line strengths of optical, infra- 
red, and ultraviolet spectral lines at the photosphere of SN 
1 987A. Concerning the earliest phase of development, the only 
calculated the optical depths of lines for temperatures above 
6000 K. He showed that the UV lines are stronger than the 
optical /IR lines, and the line strengths increase with decreas- 
ing temperature. Fransson ( 1994) obtained similar results and 
showed that even the low abundance elements in the envelope 
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and in the core have very large optical depths in the UV. In 
an atmosphere with a velocity range of 2000-20,000 kms -1 , 
typical of supernova expansion, each line occupies a wave- 
length band much greater than in a static atmosphere. The U V 
photons will be scattered from line to line by the gas in the core 
and envelope, while being redshifted each time. Some photons 
will stop resonance scattering and finally escape from the en- 
velope by finding a wavelength “window” where the scattering 
lines are neither numerous nor thick enough to scatter the pho- 
tons (Fransson 1994; Xu Sc McCray 1991 ), while some other 
photons may be destroyed by photoionization, the photoioni- 
zation of the hydrogen atoms in the Balmer continuum ( Xu & 
McCray 1 99 1 ; Xu et al. 1 992 ) . The effective opacity of such an 
expanding atmosphere, the expansion opacity, has the net 
effect of line blanketing, that is, producing a quasi-continuum 
formed by many overlapping lines (Karp et al. 1977). In the 
UV region where there are a great number of densely spaced 
resonance lines, the line opacity is spread out in the observer 
frame frequency, thus enhancing the flux reduction. Wagoner, 
Perez, & Vasu ( 1991 ) studied the expansion opacity of a hy- 
drogen-rich atmosphere similar to that of the outer layers of 
a SN II. Despite the fact that SN 1987A has a much smaller 
hydrogen envelope than other SNe II, we expect the model at- 
mosphere calculations by Wagoner et al. ( 1991 ) will also be 
valid for the SN 1987 A atmosphere because hydrogen is not 
the dominant source of opacity in the optical /UV region. They 
predict a substantial flux deficiency in the U V and that the UV 
flux deficiency increases with decreasing temperature. The ex- 
pansion opacities in the wavelength region 1000-4000 A in- 
crease by a factor of more than 100 as the temperature of the 
atmosphere drops from 1 2,000 to 5000 K. while the expansion 
opacities of the optical and IR region are not affected much by 
the drop in temperature. This suggests qualitatively why the 
UV spectrum remains opaque long after the optical spectrum 
has become optically thin. Thus at early times, if the tempera- 
ture is high enough, we should observe no UV flux deficiency 
from line blanketing. For the late-time SN 1 987A atmosphere, 
most of the lines depart significantly from LTE and the opacity 
calculations become very model-dependent. 

The development of the UV spectra of SN 1987 A in the ear- 
liest phases can be explained by the decrease in temperature of 
the ejecta and the increasing effects of line blanketing 
(Cassatella et al. 1987; Fransson et al. 1987; Kirshner et al. 
1987b; Lucy 1987). At day 1.7 (Fig. 2), we were observing 
the cooling tail of a much hotter photosphere heated by the 
interaction of the explosion shockwave with the stellar surface. 
Hydrodynamical models of the shock breakout in SN 1987 A 
by Ensman & Burrows ( 1992) showed that the UV burst oc- 
curred about 2 hr after the core collapse of the progenitor. In 
these models, the luminosity rises more than five orders of 
magnitude over an interval of a few hundred seconds and the 
effective temperature of the emitted radiation reaches about 
5 X 10 4 K. And the observed UV spectrum showed a contin- 
uum with broad features. After the first IUE observation, the 
temperature of the supernova declined rapidly as the super- 
nova atmosphere expanded and cooled. The net effect in the 
UV of the hundreds of blended iron group lines is that the UV 
continuum below 2600 A disappeared below detection by 
about 3 days after outburst. Even the lines strong enough to 
give rise to P Cygni profiles were modulated by line blending 
due to the weak lines. The observed peaks in the spectra were 


the consequences of minima, or “windows,” in selective opac- 
ity. The synthetic UV spectrum computed by Lucy ( 1987) 
confirmed that the trough observed in the IUE spectrum at day 
3 in the region 2200-2600 A is caused by the blending of more 
than 150 closely spaced lines of Fe n and Fe ni. These lines 
would be individually strong enough to be observable if not 
blended together by differential expansion. This same process 
is responsible for the spectral and bolometric evolution of clas- 
sical novae ( Hauschildt et al. 1994). 

After the dramatic changes in the early days in both the flux 
level and the spectral shape, the UV spectra showed only minor 
and slow variations at epochs after — day 1 50. We know that 
falling temperature in the early SN 1987 A atmosphere due to 
expansion will lead to increased UV opacity. However, since 
the UV region is formed by lines with very large opacities, the 
spectral shape did not respond to the changes in opacity as 
quickly as it would if the opacity were lower (Jeffery & Branch 
1990). This explains the slow evolution of the UV spectrum 
of SN 1987 A after the initial phase. While the optical spectra 
obviously became nebular (optically thin) by ~~ day 200, the 
UV spectra have remained photospheric (optically thick) due 
to the heavy line blanketing effects during the whole epoch of 
IUE observations. 

McCray, Shull, & Sutherland ( 1987) predicted that after the 
initial drop in flux, there should be an “ultraviolet renaissance” 
in the UV region for SN 1 987 A. They suggested that as the 
supernova envelope expanded and became thin, the inner de- 
bris of the explosion would be revealed, releasing high-energy 
radiation originating from a hot photosphere caused by either 
the decay of radionucleides or a newly born pulsar. Fast elec- 
trons would be produced by Compton scattering of these high- 
energy photons, and the direct excitation of H, He, and heavy 
elements in the expanding debris by these electrons would pro- 
duce radiation in UV and X-ray. The UV photons will be 
trapped by the substantial expansion opacity and come out as 
a quasi-continuum (Karp et al. 1977). They predicted a sub- 
stantial UV continuum emerging from the supernova enve- 
lope about a few months after the explosion, and the UV spec- 
trum to break up into an emission-line spectrum about one 
year later. This effect is similar to the one seen in classical no- 
vae where an increasing UV flux is observed as soon as the shell 
becomes transparent and the hot nova remnant ionizes the 
ejecta (Williams 1991; Saizar et al. 1992). However, this re- 
surgence of UV radiation was not observed by either IUE or 
HST. One explanation is that the newly formed dust mole- 
cules in the supernova ejecta may have large ultraviolet opacity 
to absorb the UV radiation and re-radiate this luminosity in 
the IR continuum, as observed in some classical novae (Gehrz 
1988 ). However, this is probably not applicable for SN 1987 A. 
As we will see in § 5.2, instead we observe a rise in the UV flux 
relative to the bolometric luminosity, even after the epoch of 
dust formation. The failed prediction of a UV renaissance 
probably lies in the fact that the UV radiation originates from 
the Doppler blanketing of emission lines from mostly iron 
group elements and not from the scattering from high-energy 
photons. 

3.2. Time Evolution 

The three-dimensional plots shown in Figure 3 (SWP) and 
Figure 4 ( LWP) illustrate vividly the time evolution of the SN 
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Fig. 3. — Three-dimensional plot of/17: SWP spectra { 1250-1975 A) ofSN 1987A in time. The combined spectrum is binned at 4 A in wavelength 
scale and 5 days in timescale. The plot is smoothed ( boxcar = 5 ) in both the time and wavelength axes. The time evolution of broad band emission feature 
A ( 1980 A) is shown. 
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Fig. 4.— Three-dimensional plot of WE LWP spectra (2000-3300 A) of SN I987A in time. The combined spectrum is binned at 6 A in wavelength 
scale and 5 days in timescale. The plot is smoothed the same way as in Fig. 3. Notice that the flux scale is larger than that of Fig. 3 by a factor of 10. The time 
evolution of broad-band emission features B ( 2100 A ). C ( 2200 A). D ( 2650 A ), E ( 3000 AT F( 3 100 A), and G ( 3247 A) are shown. 
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1987 A UV spectrum. In both figures, spectra taken before day 
8 are excluded. The vertical scales of the two figures are differ- 
ent, with that of Figure 4 larger by a factor of 10. The change 
of the shape, magnitude, and strength of various spectral fea- 
tures can be studied despite the slow rate of evolution. Similar 
display of the time evolution of the SN 1 987A spectrum is pre- 
sented in Hanuschik & Thimm ( 1990) on the evolution of the 
fine spectroscopic structures of the spectral lines obtained at 
the Bochum telescope. Different time behaviors for different 
spectral features are demonstrated in Figure 5. The evolution 
of the SN 1987A spectrum in the range 2550-3300 A, where 
most of the UV spectral features are located, is plotted. Evolu- 
tion of the strength of individual features can be studied and 
compared. For example, the evolution of the strengths of the 
two emission features at 3000 A (feature E in Figs. 4 and 5) 
and 3247 A (feature G) are similar from around day 150 to 
day 800, but the maximum in flux for the 3247 A feature at 
around day 85 is not observed in the 3000 A feature. Thus it 
seems that the maximum observed at day 85 for the 3247 A 
feature probably originates from a different ion, or from a 
different layer in the photosphere. 

Since the UV spectrum consists of broad features formed by 
the blending of many weaker lines of similar wavelengths, it is 
a difficult task to identify the individual lines. Cassatella et al. 
( 1987) identified the features between 2800 and 2950 A during 
the first 30 days of observation as contributions from the Fe 
n multiplets. Also, the observed absorption blends near 3270, 
3160, and 3015 A are identified as the resonance lines of Ti II. 
Lyubimkov (1990) studied the spectra from the Soviet 
ASTRON satellite ( to be discussed in §4.1 > nd also concluded 
that these absorption blends result from sti ng Ti n multiplets 
such as A3372 (No. 1 in the Martin, Fuhr, & Wiese 1988 list of 


Ti II multiplets), A3331 (No. 7), A3238 (No. 2), and A3080 
(No. 5). However, most other features seen in the UV spec- 
trum of SN 1987A cannot be identified with a given line or 
multiplet of any particular species. Despite the lack of specific 
identification of many features in the spectra, two trends of 
their evolutions can be deduced: first, the absorption lines be- 
come narrower in time, the flux minima around 2500, 2700, 
and 3200 A; and second, the emission maxima shift towards 
longer wavelengths, the flux maxima around 2650, 3000, 3100, 
and 3247 A. Similar trends were also observed in the optical 
spectra (Catchpole et al. 1987; Phillips et al. 1988, 1990). 
There are many factors affecting the radiative transfer in the 
late-time supernova envelope — changes in temperature and 
non-LTE conditions, energy of gamma ray photons from ra- 
dioactive nuclei, and dust formation, for example. However, 
the essential effect is probably the declining density as the in- 
verse third power of time (p — t“ 3 ) as the debris expands. Sim- 
ilar to other P Cygni absorptions, this expansion will cause the 
photosphere (above which the lines are formed) and the re- 
gions of large line opacity (the line- forming region) to recede 
to the inner, denser, more slowly expanding regions of the 
ejecta. The emission maxima will be shifted to longer wave- 
lengths because their red sides are absorbed by strong lines 
which become narrower in time. 

In order to study the redshift of the features with time, we 
have determined the peak wavelengths of six spectral features: 
1980 (feature A in Fig. 3), 2100 (feature B in Figs. 4 and 5), 
2200 (feature C), 2650 (feature D), 3100 (feature F), and 
3247 A (feature G). These features were chosen because their 
maxima are easily identifiable. A tenth-order polynomial was 
fitted to the part of the spectrum where the feature is located, 
and then we determined the wavelength maximum of the fitted 
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Fig. 5. — Same as Fig. 4 for the wavelength range 2550-3300 A is shown for clearer display of the time evolution of the emission features D-G 
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polynomial. The evolutions of the peak wavelength for the six 
features are shown in Figure 6. There are fewer data points for 
the 2100 and 2200 A features because we only included the 
spectra with good signal-to-noise in that wavelength region. 
The accuracy of the measurement is not significantly ham- 
pered by the resolution of the spectra because the supernova 
features remained quite broad. We found that except for the 
3247 A feature, the redshifting of the features from about day 
120 to day 800 appears linear with time. Second, for the 3100 
and 3247 A features, additional phases of rapid redshifting are 
observed before day 1 20. At the same time, both these features 
showed a flux maximum around day 85, which were not seen 
for the 2100, 2200, and 2650 A features ( Fig. 4). We conclude 
that the UV radiation of these features from the two epochs 
(between day 10-120, and day 120-800) probably originated 
from species in different regions of the supernova atmosphere 
which have different velocities. The 1980 A feature did not 
have maximum magnitude at day 85, yet there seems to be a 
rapid redshifting component from day 1 today 10. However, 
we suspect that the two branches of the peak wavelength evo- 
lution are the results of two unrelated blends, instead of the 
same species, because the feature is not observed between day 
10 and 55. 

The precise radial velocities of these features cannot be de- 
termined without knowing the correct rest wavelengths. How- 
ever, arbitrarily choosing an intermediate observed value as the 
rest wavelength introduces less than a 1% error in the rate of 
change of radial velocity. Velocity increases of 3. 74, 4.09, 3.95, 
3.85, and 4.38 km s 1 day " 1 are determined from linear fits to 
the emission maxima of the 1 980, 2100, 2200, 2650, and 3 100 
A features respectively from day 120 to 800. The average rate 
of velocity increase is found to be 4.00 km s 1 day 1 . This 
result is surprisingly similar to the results obtained by fitting 
comparable data of the Ha P Cygni absorption minimum from 
day 120 to day 526. A straight line through the optical veloci- 
ties of the Cerro Tololo Inter- American Observatory (CTIO) 
data shows that the rate of radial velocity change is 4.23 km s 1 
day 1 (Phillips etal. 1988, 1990). 

Further comparisons with optical results are provided by the 
development of the 3247 A feature. Employing 3247 A as the 
rest wavelength and assuming a radial velocity for the progeni- 
tor of +286 km s _1 (Wampler, Richichi, & Baade 1989), the 
variation of radial velocity with time can be determined. The 
radial velocity evolution of the 3247 A feature is compared to 
that of the Fe n 5018 absorption, and the results are presented 
in Figure 7. The two sets of data agree extremely well starting 
from the first day of the detection of the 3247 A feature at day 
6.0. The agreement continued until day 45, after which they 
started to differ. The result suggests that the species contribut- 
ing to the 3247 A emission blend had the same velocity as Fe II 
from day 6 to day 45, implying that they originated from the 
same layer in the ejecta. There is also the possibility that the 
3247 A blend in fact has a Fe li origin from day 6 to day 45, 
while other components in the blend, possibly Ti il 
(Lyubimkov 1990), dominated afterward. 

3.3. Combined Optical-UV Spectra 

To better understand the evolution of SN 1987A, the ultra- 
violet WE observations have been combined with the optical 


spectra taken at CTIO. From day 1 to day 1 30, the optical spec- 
tra covered the wavelength range 3700-7100 A at a FWHM 
resolution of 5 A (Phillips et al. 1988). From day 198 up to 
day 805, observations with wider wavelength coverage (3200- 
11,000 A) but lower resolution (11-16 A) were obtained 
(Phillips et al. 1990). An observation on day 1050 obtained 
at slightly higher resolution ( 10 A FWHM) at the CTIO 4 m 
telescope is also included (Suntzeff et al. 1991 ). The optical 
spectra are combined with IUE data taken within 0.4 days be- 
fore day 200, and within 2 days after day 200, when there are 
fewer optical and UV spectra available. The 32 UV + optical 
combined spectra are shown in Figure 8. Note that all the 
CTIO optical spectra taken prior to day 805 contained con- 
tamination from the field stars near SN 1987 A, namely stars 2 
and 3. As the supernova faded, the percentage of contamina- 
tion from these field stars increased. The corrections for con- 
tamination are significant for spectra obtained after about day 
600. For example, at day 600 about 35% of the flux at U and 
17% at / comes from the field stars ( Phillips et al. 1990). Spec- 
trophotometry was obtained later for both star 2 and star 3 
( Suntzeff et al. 1991; Walborn et al. 1 993 ) . However, in Figure 
8, only the day 1050 spectrum has these two stars removed, 
where almost 70% of observed T/-band flux is from stars 2 and 
3. The optical spectra are corrected for interstellar extinction 
with the average extinction law of Cardelli, Clayton, & Mathis 
( 1989). To be consistent with the correction of WE data, we 
use E(B - F) lotal of 0.15 and R v of 3.1. Since the UV flux 
decreased at an enormous rate in the first few days of observa- 
tion, a small time difference can lead to a big difference in ob- 
served flux. For example, in the earliest spectrum (day 1.7), 
the UV spectrum was collected 0.2 day earlier than the CTIO 
one. During the time the rate of change of U V flux in the near 
UV region of 3000-3300 A is AFJ At ~ -0.1 dex/day. Thus 
the U V flux level should be lowered by —0.02 dex for compar- 
ison with the optical data. However, this effect is only signifi- 
cant in the very early days of observations (before ~~ day 6) 
when the UV flux was plummeting, with A FJ At as large 
as — -0.6 dex/day in the overlap region (as discussed further 
in § 5 ). 

In Figure 9, A/* is plotted with A ( 1 150-1 1,000 A) for the 
combined U V T optical spectra taken at selected dates to illus- 
trate the energy budget in the UV and optical regions for SN 
1987 A. On the first day of observation at day 1.7, the energy 
output of the supernova is more or less evenly distributed from 
about 1 300 A to the red end of the optical spectrum. (The steep 
drop blueward of 1300 A is not an instrumental artifact.) The 
radiation energy in the UV range 1250-3300 A, £ uv , is 38% of 
that in the optical + UV range 1250-7150 A, Euv+opt.cai ■ This 
is compared to the value of 35% if the energy were evenly dis- 
tributed in all wavelengths. Within the next few days, the dis- 
tribution of the radiation energy changed swiftly. By day 4.6, 
most of the energy radiated from the supernova is from the 
optical regime and the value of ifuv/fuv+optjcai dropped to 
2.4%. The UV flux level continued to decrease while the optical 
flux was increasing, releasing the hydrogen recombination en- 
ergy. At day 44, the energy level coming out from the optical 
region is about 10 4 times of that in the UV and the ratio £uv/ 
£uv+o P ticai dropped to 0.05%. After reaching the maximum at 
day 85, the optical flux started to decay linearly, while the UV 
remained steady (§ 5.2, Fig. 16). At day 198, the optical flux 
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Days after outburst 

Fig. 7. — Comparison of the measured radial velocity of feature G 
(3247 A) with that of the Fe u 5018 absorption line measured at CTIO 
( Phillips et al. (988). The radial velocity of the progenitor is assumed to be 
+286 km s~ l . 


was lower that seen on day 44. But at the same time, the UV 
flux was higher than that of the level at day 44. Therefore, while 
the energy emitted at optical and UV wavelengths differed by 
about one order of magnitude at day 423 (£ uv /£uv+opticai — 
2%), they are nearly at the same level at day 715 (E uv / 
^iiv+opticai — 14%). As the supernova aged, relatively more en- 
ergy of SN 1987A comes from the ultraviolet region compared 
to the optical. 

4. COMPARISONS TO OTHER ULTRAVIOLET 
OBSERVATIONS 

4.1. ASTRON Satellite 

The UV spectra from IUE are compared with the results 
from the Soviet ASTRON satellite. The “Astrophysical Station 
ASTRON” was launched on 1983 March 23, and it continued 
to function until mid- 1988. The detailed characteristics of the 
satellite were described by Boyarchuk et al. ( 1 986 ) . It consisted 
of two instruments, one of which is a f/ 10 Ritchey-Chretien 
ultraviolet telescope. The detector consists of a three-channel 
photoelectric scanning spectrometer covering the spectral re- 
gions 1050-3500 A. 

Fifteen ultraviolet spectra were taken with ASTRON on SN 
1987A covering the spectral range 2600-3500 A from day 9 to 
day 394 (Boyarchuk et al. 1987; Lyubimkov 1990). With the 
supernova observed through an entry aperture 1 ' in diameter 
at a resolution of 28 A, the ASTRON spectra should contain 
flux from star 2 and star 3 in the field of SN 1 987 A. However, 
in the spectral region where IUE and ^ 577? CW overlap ( 2800- 
3300 A, low signal-to-noise for the 2600-2800 A data ) , the flux 
level of the background stars, 4.4 X 10~ 14 ergscnT 2 s~ l A"\ is 
less than 4% of the flux of the supernova during the ASTRON 
observation epochs. Thus the background stars did not have 
significant flux contribution to the ASTRON spectra. These 
spectra and the corresponding IUE data as observed (before 
reddening corrections) at similar epochs are presented in Fig- 
ure 10. Each successive ASTRON spectrum has been shifted 
down 0.2 dex relative to the preceding one to avoid overlaps. 


All IUE spectra displayed in Figure 10 are smoothed to the 
resolution of the ASTRON observations. Each successive IUE 
spectrum is also shifted down by 0.2 dex, except for day 1 20.5, 
which is downshifted 0.4 dex relative to the preceeding spec- 
trum. The two sets of observations generally agree well on the 
spectral shapes and the relative strengths of the spectral fea- 
tures in the overlapping region. However, the absolute flux 
measured for these features were systematically different, with 
the IUE results higher by ~~0.02-0. 1 dex. The difference prob- 
ably results from the different UV calibration techniques for 
the two satellites. 

4.2. VBLUW Photometry 

The IUE results can also be compared with the ground- 
based VBLUW Walraven photometry. Walraven photometry 
for SN I987A was measured at the 9 1 cm Dutch Telescope of 
the European Southern Observatory (ESO) at La Silla, Chile. 
The fF-band of the system, with effective wavelength of 3236 
A and FWHM bandwidth of 1 57 A, provides a useful compar- 
ison to the near-ultraviolet results of IUE. Photometry results 
from day 1 to day 393 ( 1988 March 22) were obtained, and 
some of the early results have been discussed by Cristiani et al. 

( 1987) and Pel et al. ( 1987), The IUE spectra, not corrected 
for reddening, are convolved with the transmission function of 
the fF-band as provided by J. W. Pel. The fF- band photome- 
try can thus be directly compared with the mean flux density 
F X (P), which is given by the expression 

f aC J P(X)/x(X)XJX 

= , (l) 

P(\)\d\ 

Jo 

where P(\) is the filter passband andjx(X) is the spectroscopic 
flux (Horne 1988). Equal weight is given to equal photon 
number rather than equal photon energy in equation ( 1 ). In 
Figure 1 1 , both the IUE synthesized mean flux results and the 
fF-band photometry data are presented. This provides a con- 
venient link between ground-based photometry and the IUE 
data. The accuracy of the B'-band photometric data on SN 
1 987A is about 0.005 dex ( Pel et al. 1 987 ) . The two sets of data 
agree extremely well — the rms difference between the two data 
sets is 0.02 dex. We conclude that the absolute calibration of 
the LWP spectra at ~3200 A is correct within an error of ^5% 
and the calibration procedures used by the ASTRON generate 
results in which the mean flux appears lower by ~~ 15%. 

5. THE UV LIGHT CURVE OF SN 1987A 

5.1. UV Photometry 

The IUE spectra have been integrated over a limited number 
of broad-wavelength intervals so that time variations of the ob- 
served flux can be studied. The detection range of the LWP is 
divided into three wide bands: short, S: 1975-2500 A; middle, 
M: 2500-3000 A; long, L: 3000-3300 A (not to be confused 
with the infrared M and L bands). The flux collected from the 
SWP camera ( 1 1 50- 1 975 A ) is also integrated. We did not cor- 
rect the SWP photometry for the circum stellar lines because of 
their small contribution to the total flux (<2%). The wave- 
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Fig. 8.— (a )-(//) The combined UV and optical spectra ofSN 1987 A. Notice that both the horizontal axis (wavelength) and the vertical axis (flux) are logarithmic. The optical spectra are obtained 
from CTIO ( Phillips et al. 1988. 1990; SuntzefFet al. 1991 ). The optical spectra are corrected for interstellar extinction with the Cardelli et al. ( 1989) law with E(B ~ V ) loU) = 0. 15 and R y = 3.1. 
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Fig. 9. — Plots of X/ x vs. X ( 1 150-1 1,000 A) at selected dates to show the energy distribution in the spectra of SN 1987A 


length band from 1 199 to 1227 A is excluded from the SWP 
integration because of strong geocoronal Lya emission. The 
UV light curves from these wavelength bands for the first 1 5 
days after the explosion are presented in Figure 12. The inte- 
grated mean flux density of the waveband from X i to Xyris de- 
termined by the formula 


Fx = 


f / w x A(X)X^X 

A 

fV 

w x \d\ 


( 2 ) 


where w x is the Lya notch filter such that w x = 0 from 1 199 
to 1227 A and w x = 1 elsewhere. The flux from SN 1987 A 
plummeted over the first few days of observations for all the 
UV wavelength bands. The decrease continued throughout the 
1 5 day period, though it was much slower after day 8. The time 
development of the integrated flux from the SWP band is very 
similar to that of the S band, while the M and L bands showed 
different temporal developments. During the initial drop, the 
short-wavelength bands declined much more rapidly than the 
long-wavelength bands. During the first 5 days of IUE obser- 
vations, the integrated flux of the wavelength bands shortward 


of 2500 A decreased by a factor of 2000, while the flux from 
the wavelength bands longward of 2500 A decreased by a factor 
of 100. This steep drop in UV flux is accompanied by a large 
variation of the UV spectral characteristics as discussed in § 3. 

The dramatic initial drop in U V flux can be explained by the 
rapid temperature decrease in the photosphere of the super- 
nova, and also the corresponding increase in expansion opacity 
with the falling temperature (§ 3.1). We are probably seeing 
the cooling tail of a much hotter photosphere which was pres- 
ent on the day of outburst when the shock emerged from the 
stellar interior of Sk -69 202. The effects of the falling photo- 
spheric temperature are illustrated in Figure 1 3 where the com- 
bined optical + UV spectra of the SN 1987 A for day 1.7 and 
day 4.6 are fitted by blackbody spectra of temperatures 1 1 ,700 
and 7500 K, respectively. The measured temperatures are de- 
termined by fitting the IUE and optical /IR colors (Hamuy et 
al. 1988). The scale of the blackbody spectrum was determined 
by demanding the integrated fluxes of the observed and black- 
body spectra to be equal for the range 5000-6200 A 
(approximately the Kband). For the day 1 .7 spectrum, all the 
flux points from UV to the M band are used for the fit. After 
that, the UV and U-band fluxes are systematically lower than 
a blackbody and thus only the B to L magnitudes are used for 
the fit to the day 4.6 spectrum. The initial temperature drop 
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Wavelength (Angstrom) 

Fig. 10.—' The UV spectra taken by the ASTRON satellite compared with the IUE data at similar dates. The IUE data is smoothed to the resolution of 
the ASTROS’ fax*. ( 28 A). Each successive spectrum is shifted downwards by 0.2 dex for clarity of display. The only exception is the IUE spectrum taken at 
day 1 20.5, which is shifted down by 0.4 dex. Neither set of data is corrected for interstellar extinction. 



Fig . 1 1. — Comparison of the Walraven IT -band photometry (A** = 
3236 A; FWHM = 157 A) with the convolved fT-band IUE spectropho- 
tometry results. 


has two effects on the evolution of the UV flux. First, decreas- 
ing temperature causes the blackbody emission maximum to 
shift fromthe UV regime ( ^2500 A) to the blue optical region 
(~~3900 A). Second, the falling temperature vastly increases 
the expansion opacity of the supernova atmosphere ( Wagoner 
et at. 1991 ). They showed that at 5000 K, the UV expansion 
opacity is about 10 4 times that of the optical/IR. At day 1.7, 
the UV flux level slightly exceeded that of the blackbody fit. 
However, in two days, the situation quickly reversed and SN 
1987A showed definite UV flux deficiency relative to black- 
body fits to the optical. 

We tried to estimate the effects of increased opacity in UV 
during the earliest stage of the supernova. The ratio between 
the integrated UV flux of the observed IUE spectrum and the 
fitted blackbody spectrum is computed for the first 1 5 days of 
observations. The temperature fits from Hamuy et al. ( 1988) 
are used to construct the blackbody spectra. Because of the 
similar time development of the S and SWP band flux in the 
early phase (Fig. 12), they are combined for the calculation 
(the SWPTS band). The results for the integrated flux ratio of 
the SWP+S, M, and L bands, along with the full IUE band. 
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Fig. 12. — UV light curves of SN 1987A of various wavelength bands 
for the first 1 5 days after outburst. 


are presented in Figure 14. We notice that the UV deficiencies 
increased with time for all the wavebands. While the UV flux 
level is about 40% above the blackbody radiation at day 1 .7, it 
dropped to be only 10% of the blackbody flux around day 5, 
and 1% around day 14. Second, the UV deficiency is more sig- 
nificant in the shorter wavelength band. The SWP+S band is 
more underluminous to a blackbody than the L band by about 
one order of magnitude between day 4 and 1 4. This is probably 
caused by the enormous line-blanketing effects in the wave- 
length region below 2600 A with the presence of many strong 
Fe II lines. 

Figure 12 shows that the UV flux decreases more rapidly at 
shorter wavelength during the early days of observations. This 
can also be seen in Figure 1 5, where the rates of decrease of the 
IUE flux, - A FJ At (expressed in unit of magnitude day - 1 ), in 
various wavebands are plotted against time. The earliest UV 
flux decay rates measured at day 2 are 1 .93, 0.56, and 0.20 mag 
day " 1 , for the SWP+S, M, and L wavebands, respectively. The 
rates of change of the optical (7, B y and V bands photometry 
are also included in the figure. No interstellar extinction cor- 
rection has been applied to the IUE data so that direct compar- 
ison between the rates of flux change in the optical and in the 
UV is possible. From day 3.8 through day 813, the 
UBV(RI) kc photoelectric observations taken at CTIO 
(Hamuy & Suntzeff 1990) are used, while prior to day 3.8, the 
UBVRI photoelectric observations by Shelton ( 1993 ) are used 
to compute the rate of change of flux in £7, B, and V. Discrep- 
ancies exist between the two data sets, mostly because of the 



Wavelength (A) 

Fig. 13. — Comparingthe earliest UV + optical spectrum of SN 1987A 
at day 1 .7 (extinction corrected) with the blackbody spectrum of tempera- 
ture 1 1 700 K. The BB temperature is obtained by fitting the optical and IR 
colors ( Hamuy el al. 1988). The spectrum at day 4.6 and the blackbody 
spectrum of 7500 K is also presented. 



Fig . 1 4.— -Ratios between observed flux and the fitted blackbody spec- 
trum in various UV wavelength bands. 
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Fig. 15. — Rate of change of flux in time for various UV wavelength 
bands and optical photometry. The optical photometry data are from Shel- 
ton ( 1993) and the CTIO observatory (Hamuy Sc Suntzeff 1990). The 
smooth lines show the rates of change of U V flux for the fitted blackbody 
spectrum with temperature T(t)= T 0 [l + k exp (-at)], where T 0 = 5900 
K, A; = 3.0 and a = 0.54. The SWP+S {solid line), M ( dotted line), L 
{broken line) results are shown. 


different filter passbands used. However, for the first few days 
after outburst, we estimate that the error involved should be 
less than —0.03 mag. The trend of slower decay rates with 
longer wavelength extends to the U and B bands in the optical 
regime. In fact, the Kband flux, along with R and I band fluxes 
(not plotted here), increases during the early days of observa- 
tion. The decay rate for the SWP+S band reached a maximum 
of 2. 1 3 mag/day at day 3. 1 , while the decay rates for both the 
M and L bands reached maxima 1 .3 days later at day 4.4, with 
the value of 1.59, and 1.51 mag/day, respectively. After that, 
the decline in UV flux became slower with time in all wave- 
bands. Despite the big differences in the initial decay rates ob- 
served in the various wavelength bands, all the U V bands decay 
at similar rates starting from around day 6.5, reaching a value 
of —0. 1 5 mag/ day at day 1 4. The slow down of flux decay in 
all UV bands may probably be explained by the diffusion of 
energy which was deposited by the passage of the initial shock. 
However, the exact understanding of the phenomenon re- 
quires detailed modeling of the UV radiation from the atmo- 
sphere at the early stages. 

The evolution of — A FJ A/, for a fitted blackbody is also pre- 
sented in Figure !5. We estimate AFJAt by dF x /dt = (dFJ 
dT)(dT/dt). SFJdT is determined by differentiating the 
Planck’s function, while dT/dt is calculated by fitting the early 
temperature data from CTIO ( Hamuy et al. 1 988 ) to a simple 
analytic expression T{t) = T 0 [l + k exp(-aO). The values 
T 0 = 5900 K, k - 3.0 and a = 0.54 are used for the fit. This 
model is able to reproduce the general temporal behaviors of 
-AFJAt observed in all UV wavebands. The peaks of flux 
change rate occur earlier for the shorter wavelength bands ( day 
1.4 for SWP+S band) than the long-wavelength ones (day 1.8 
for M and L bands). However, the peaks generated by the 
model occurred about 2 days earlier than the observed values 
and the model peaks for the M and L bands are too broad com- 
pared to the data. This calculation reaffirms the notion that the 


initial drop in temperature alone is not sufficient to completely 
explain the early behavior of the UV flux. 

The light curves from day 1 to day 1 248 for all the UV wave- 
length bands are illustrated in Figure 16. The integrated IUE 
mean flux densities from equation ( 2 ) , without reddening cor- 
rections, are transformed to magnitude with the formula 

magnitude = -2.5 log 10 ( integrated flux ) + constant , ( 3 ) 

where the value of constant is chosen for clarity of display. Fol- 
lowing the fast decay of UV flux in the early days of observa- 
tion, each of the UV light curves ceased to decrease when a 
minimum at day 30 was reached. For the two UV wavelength 
bands longward of 2500 A (M and L), observed flux reached a 
maximum at —day 85 before reaching minimum again at 
—day 125. The UV flux from SN 1987A then increased slowly 
to another maximum at —day 350 before its ultimate decline. 
The rate of decay increased from —day 350 to —day 650, and 
remained constant until —day 900. During this constant decay 
epoch, the flux decline is faster for longer wavelengths ( 0.0109 
and 0.0129 mag/day for M and L band respectively). It is 
difficult to determine the rate of decay of the light curves after 
day 900 as the supernova became dimmer. But the M-band 
light curve, which has the highest signal to noise at late times, 
suggests that the UV light curves decline at a slower rate from 
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Fig. 16. — Light curves of SN 1 987A for both the UV and optical wave- 
length bands from day 1 to day 1 300. The optical photometry data are 
from Shelton ( 1993), Hamuy & Suntzeff ( 1990), and Walker Sc Suntzeff 
( 1991 ). The UV broad-band spectrophotometry results are shifted verti- 
cally for clarity of display. 
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—day 950 until the final data point at day 1248. For the two 
wavelength bands shortward of 2500 A (SWP and S), the tem- 
poral development of the light curve is similar, except that the 
second maximum near day 85 is not observed. The rate of de- 
cay for the constant decay epoch from day 650 to day 800 is 
0.0092 mag/day for both wavebands. Moreover, the similarity 
between the SWP and S band flux evolution as noticed in the 
early days of observation continued in the late epoch. The two 
light curves did not start to differ in shape until —day 900, 
when additional radiation from the variable star 3 might be 
present in the SWP spectra (to be discussed later). 

In Figure 16, the UV spectrophotometry results are com- 
pared with the optical light curves from the CTIO data, which 
include the photoelectric photometry results from day 3.8 to 
8 1 3 ( Hamuy & Suntzeff 1 990 ) , and the CCD UB VRI photom- 
etry measurements from day 680 to 1469 (Walker & Suntzeff 
1991 ). Early photometric results by Shelton (1993) from day 
1 .0 to 3.7 are also included despite the slight difference between 
the filters used. We can see that with decreasing wavelength 
from I band to far-UV, the magnitude of the light curve’s sec- 
ond maximum around day 85 gradually diminishes. This max- 
imum in fact disappeared for wavelength bands below 2500 A. 
On the other hand, after day 125, the decline of V y R y and / 
light curves are within 1 5% of the e-folding time of the radio- 
active decay of 56 Co (111.3 d), and the B band had a less rapid 
decline of — 1 60 days ( Suntzeff et al. 1988). However, this de- 
cline is not observed in the U band and UV photometric mea- 
surements. Instead, the development of a third maximum 
around day 350 can be seen with decreasing wavelength. Both 
the disappearance of the second maximum of the UV light 
curves around day 85 and the appearance of the third maxi- 
mum maximum around day 350 with decreasing wavelengths 
are probably related to the changes of the ionization and exci- 
tation stages of the iron group elements. However, a plausible 
interpretation of the UV light curves is not possible without 
correct modeling of the opacities with the changing tempera- 
ture, density, and composition of the photosphere. 

5.2. UVOIR Total Bolometric Luminosity 

The total integrated UV light curve ( 1 150-3300 A) of SN 
1987A is shown in Figure 17. The same reddening correction 
as described in § 2.3 is applied to the data, i.e., E(B - K) tolal = 
0.15. Only epochs with both SWP and LWP spectra taken 
within 5 hr are investigated. The spectra are combined for cal- 
culations and a total of 96 measurements are obtained. Also 
shown in Figure 17 is the total ultraviolet-optical-infrared 
( UVOIR) bolometric light curve of SN 1987 A. The CTIO and 
ESO groups have previously computed the UVOIR bolometric 
light curve based on broadband U to M photometry for the 
early days of observations (Suntzeff & Bouchet 1990), and U 
to Q 0 photometry for the late-time observations ( Suntzeff et al. 
1991; Bouchet et al. 1991 ). The ESO and CTIO results differ 
from each other after day 6 1 6 because of the differences in the 
1R flux measurements. In these compilations, the UV flux was 
included in the total bolometric flux calculations for only the 
first 14 days. After day 14, the UV contribution was assumed 
to be negligible, and thus only the optical/ IR luminosities, 
Lqi R , are presented. While this omission was justified before 
—day 200, the UV flux contributes an increasing fraction of 



Fig. 17. — Revised UVOIR bolometric light curve of SN 1987A with 
UV data included from day I to 1352. Results from both CTIO ( filled 
circles) and ESO (open circles) are shown. Errors are less than the symbol 
size if not shown. The total UV light curve measured from IUE data 
(squares) and the spline fit to the UV light curve used for calculations of 
the bolometric luminosity are plotted. Extrapolations of the spline fit be- 
yond the final UV flux points are also shown (crosses). 


the total bolometric luminosity as time progresses. Therefore, 
we tried to add the UV flux contributions after day 14 to the 
bolometric data after day 14 from both the CTIO and ESO 
groups. Moreover, the IUE data that had been included for 
bolometric calculations before day 14 were not processed for 
the corrections described in § 2 such as appropriate combina- 
tion, background subtraction, camera degradation, and so on. 
Therefore, the early bolometric fluxes before day 14 are recal- 
culated with our new processed set of IUE data. The revised 
bolometric luminosity data, including the contributions of U V 
at all epochs, for SN 1987A are given in Table 4A and 4B, 
for the results from the CTIO and ESO groups, respectively. A 
distance modulus to the LMC of 1 8.5 is assumed. 

Before day 1 4, we started with the UV to Af-band bolometric 
luminosity data as described in Suntzeff & Bouchet ( 1990). 
The separate UV and optical/ IR components of the bolomet- 
ric flux used in the calculation are obtained (N. Suntzeff, pri- 
vate communications). We fit a spline to our new UV light 
curve to estimate the UV fluxes on days with OIR measure- 
ments. The interpolated UV flux is added to the measured OIR 
flux to obtain the UVOIR bolometric luminosity. The adjust- 
ment is —0.06 dex for the first day of observation (day 1.14), 
and it decreases with time. By day 10, the adjustment to the 
bolometric luminosity became <0.001 dex. Notice that the 
first two OIR measurements (day 1.14 and 1.51 ) were made 
before the earliest IUE observation (day 1 .66). Thus extrapo- 
lation of UV data is required for the bolometric flux calcula- 
tions for the two dates day 1.14 and 1.51. We estimated that 
the extrapolation will not induce large enough error to change 
the error quoted for the bolometric luminosity (0.02 dex). 

After day 14, we started with the U to Qo bolometric lumi- 
nosity data of the CTIO (Suntzeff et al. 1991) and ESO 
(Bouchet et al. 1991 ) groups. Again, we interpolate the UV 
light curve with a spline to estimate the UV fluxes on days with 
optical /IR measurements. This method works fine for epochs 
before —day 900. However, after day 900, as the supernova 




No. 1, 1995 


WE UV OBSERVATIONS OF SN I987A 


253 


TABLE 4 


A. Bolometric Luminosity of SN 1987A from CTIO (Suntzeff& Bouchet 1990; Suntzeff et al. 199I)and 
the New UVOIR Bolometric Luminosity with WE Results Included 


Day 

CTIO 

l°g (Luvoir)* 
(ergs s“‘) 

log(Luv) 
(ergs s ') 

New 

log (Luvoir) 
(ergss -1 ) 

a log (Lboi) 

Day 

CTIO 

log (Luvoir)* 
(ergs s"‘) 

log (L uv ) 
(ergs s’ 1 ) 

New 

log (Luvoir) 
(ergss -1 ) 

(HogCL**) 

1.14 

41.687 

41.446 

41.628 

<0.02 

103.70 

41.829 

38.514 

41.829 

<0.02 

1.51 

41.582 

41.275 

41.546 

<0.02 

104.70 

41.812 

38.500 

41.812 

<0.02 

1.85 

41.493 

41.118 

41.476 

<0.02 

106.70 

41.787 

38.473 

41.787 

<0.02 

2.70 

41.373 

40.722 

41.372 

<0.02 

108.70 

41.755 

38.446 

41.755 

<0.02 

3.65 

41.294 

40.226 

41.288 

<0.02 

1 10.70 

41.722 

38.421 

41.722 

<0.02 

4.70 

41.245 

39.613 

41.237 

<0.02 

111.70 

41.706 

38.409 

41.706 

<0.02 

5.64 

41.219 

39.288 

41.215 

<0.02 

115.70 

41.630 

38.367 

41.630 

<0.02 

6.82 

41.21 1 

38.972 

41.207 

<0.02 

122.70 

41.544 

38.316 

41.544 

<0.02 

7.81 

41.206 

38.806 

41.203 

<0.02 

132.10 

41.501 

38.311 

41.501 

<0.02 

8.79 

41.215 

38.679 

41.213 

<0.02 

157.10 

41.390 

38.358 

41.390 

<0.02 

9.80 

41.225 

38.573 

41.225 

<0.02 

158.10 

41.389 

38.362 

41.389 

<0.02 

10.79 

41.237 

38.491 

41.237 

<0.02 

159.10 

41.381 

38.366 

41.381 

<0.02 

13.80 

41.293 

38.348 

41.293 

<0.02 

164.10 

41.361 

38.383 

41.361 

<0.02 

14.80 

41.309 

38.326 

41.309 

<0.02 

195.00 

41.230 

38.448 

41.231 

<0.02 

19.70 

41.371 

38.261 

41.371 

<0.02 

198.00 

41.214 

38.452 

41.215 

<0.02 

24.80 

41.438 

38.220 

41.438 

<0.02 

201.00 

41.202 

38.456 

41.203 

<0.02 

25.80 

41.449 

38.229 

41.449 

<0.02 

260.80 

40.950 

38.504 

40.952 

<0.02 

30.80 

41.504 

38.234 

41.504 

<0.02 

273.90 

40.894 

38.508 

40.896 

<0.02 

32.80 

41.525 

38.219 

41.525 

<0.02 

280.90 

40.869 

38.511 

40.871 

<0.02 

33.80 

41.536 

38.214 

41.536 

<0.02 

293.00 

40.812 

38.517 

40.814 

<0.02 

34.80 

41.544 

38.210 

41.544 

<0.02 

308.80 

40.744 

38.521 

40.747 

<0.02 

35.80 

41.556 

38.208 

41.556 

<0.02 

316.80 

40.712 

38.519 

40.715 

<0.02 

36.80 

41.569 

38.208 

41.569 

<0.02 

348.80 

40.576 

38.517 

40.580 

<0.02 

37.80 

41.584 

38.209 

41.584 

<0.02 

361.80 

40.518 

38.517 

40.522 

<0.02 

39.70 

41.607 

38.215 

41.607 

<0.02 

377.80 

40.450 

38.510 

40.455 

<0.02 

42.70 

41.640 

38.231 

41.640 

<0.02 

399.70 

40.351 

38.473 

40.357 

<0.02 

43.70 

41.653 

38.238 

41.653 

<0.02 

433.00 

40.186 

38.434 

40.194 

0.02 

44.70 

41.668 

38.245 

41.668 

<0.02 

455.70 

40.100 

38.396 

40.109 

0.02 

46.70 

41.688 

38.260 

41.688 

<0,02 

468.00 

40.040 

38.385 

40.050 

0.02 

47.70 

41.699 

38.268 

41.699 

<0.02 

518.10 

39.803 

38.295 

39.816 

0.02 

58.70 

41.796 

38.382 

41.796 

<0.02 

526.10 

39.761 

38.264 

39.775 

0.02 

59.70 

41.805 

38.395 

41.805 

<0.02 

648.00 

39.25 

37.92 

39.27 

0.04 

68.70 

41.861 

38.494 

41.861 

<0.02 

670.00 

39.14 

37.84 

39.16 

0.04 

73.70 

41.880 

38.516 

41.880 

<0.02 

935.00 

38.12 

36.86 

38.14 

0.07 

74.70 

41.885 

38.458 

41.885 

<0.02 

1055.00 

37.91 

36.60 

37.93 

0.10 

79.70 

41.899 

38.553 

41.899 

<0.02 

1115.00 

37.65 

36.50 

37.68 

0.11 

80.70 

41.901 

38.584 

41.901 

<0.02 

1143.00 

37.59 

36.46 

37.62 

0.12 

98.70 

41.867 

38.569 

41.867 

<0.02 

1352.00 

37.12 

36.10 

37.16 

0.21 

101.60 

41.844 

38.540 

41.844 

<0.02 







B. Late-time Bolometric Luminosity of SN 1 987 A from ESO( Bouchet etal. 1991) and 





the New UVOIR Bolometric Luminosity with WE Results Included 




ESO 


New 



ESO 


New 



log (Lqir) 

log (Luv) 

log (Luvoir) 



log (Lqir) 

log (Luv) log (Luvoir) 


Day 

(ergs s _l ) 

(ergs s _l ) 

(ergs s’ 1 ) 

<r log (Lboi) 

Day 

(ergs s _l ) 

(ergs s *) 

(ergss *) 

* log (Lb d) 

615.91 

39.408 

38.013 

39.425 

0.007 

903.10 

38.383 

36.934 

38.398 


634.80 

39.352 

37.964 

39.369 

0.005 

991.56 

38.300 

36.689 

38.311 


667.85 

39.177 

37.850 

39.197 

0.016 

1030.50 

38.303 

36.633 

38.312 

mm 

701.85 

39.075 

37.771 

39.096 

0.018 

1054.30 

38.258 

36.598 

38.267 

■miViVM 

728.75 

38.932 

37.617 

38.953 

0.033 

1081.13 

38.080 

36.556 

38.093 

0.090 

770.71 

38.786 

37.451 

38.806 

0.049 

1 1 10.08 

37.920 

36.510 

37.937 


791.67 

38.622 

37.386 

38.646 

0.042 

1316.00 

37.680 

36.158 

37.693 



1 For time before day 14, the UVOIR bolometric luminosity Luvoir data is tabulated. After day 14, the OIR luminosity Loir is tabulated. 


dimmed, an increasing fraction of the observed UV flux is 
from the background stars, and it becomes hard to correctly 
determine the UV contribution to the UVOIR bolometric lu- 
minosity at the very late epochs ( £ day 900). 

In Figure 17, we observed a bump in the UV light curve 


between about day 1 000 to i 1 50. The bump remains even after 
careful removal of star 4 as described in § 2.2. Most of the UV 
flux we measured at late times comes from star 2 and star 3. 
Since star 3 is a Be star, it is possible that the additional flux is 
from star 3, which is known to be a variable. With the back- 
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ground flux determined by averaging spectra over the epoch 
from day 1301 to day 1 567 ( § 2.2 ) , our net U V flux measure- 
ment is vulnerable to variations in the background stars on 
shorter time scales. This idea is strengthened after studying the 
UBV(RI) kc photometry of star 3 taken by Walbom et al. 
( 1993) which is reproduced in Figure 18. Magnitudes in all 
optical bands change with time and variations of the order of 
—0.3 mag are observed. During the period day 1080-1220, 
star 3 was exceptionally bright across all wavelength bands. 
The coincidence in time of maximum optical emission from 
star 3 with the epoch of extra UV flux provides evidence that 
this extra flux is mostly from star 3. We estimate that a varia- 
tion of 15% in the UV flux from star 3 would be sufficient to 
explain the detected increase. Be stars are known to be vari- 
ables in both UV and optical regions. A 50% decrease in UV 
flux has been detected in a Be star ( Barylak & Doazan 1986), 
thus a 1 5% variation over 1 year seems plausible. 

The other potential explanation for the bump in the UV 
fight curve is that the extra flux comes from the supernova it- 
self. A possible evidence is that the increase of the UV flux 
around day 1070 is at about the same time as the increase in 
ESO bolometric luminosity observed around day 1035. At this 
very late time, the UVOIR bolometric luminosity was domi- 
nated (—80%) by infrared radiation longward of 10 *im. In 
fact, the ESO 10 ^m and 20 ^m broadband photometric fight 
curves display bumps of —0.8 mag near day 1030 ( Bouchet et 
al. 1991 ). However, the optical light curves of the CTIO data, 
for which the supernova was resolved from its neighboring 



Fig. 18. — The UBVRI photometry of Star 3 taken by Walbom et al. 
(1993). 


stars, did not show bumps at a similar epoch (Walker & Sunt- 
zeff 1991). We cannot think of any physical process originated 
from the supernova which can explain the increase of flux in 
far-IR and UV at the same time without increasing the optical 
output. Thus it seems unlikely the extra flux observed is intrin- 
sic from the supernova. 

With the ESO 10 and 20 /im broadband photometric results 
taken with a5' aperture, we suggest that the observed increase 
of flux near day 1030 may be explained by the variable IR be- 
havior of the Be stars. While studies of the far-IR variability of 
Be stars are not available, studies of the near-IR variability of 
Be stars show that the fraction of Be stars that are variables 
increases with wavelength and the amplitude of variations in- 
creases with wavelength ( Dougherty & Taylor 1 994 ) . Walbom 
et al. ( 1993) observed variations in the near-IR photometry of 
star 3 of — I mag between day 1 150 and day 1870. Thus we 
suggest it is possible that star 3 also vary in magnitude of — 1 
mag in the far-IR bands to explain the bumps observed in the 
fight curves near day 1030. 

Assuming that the extra flux in the bump of the UV fight 
curve seen between day 1000-1 150 is from the variable star 3, 
we tried to estimate the fraction of the observed UV flux that 
was originated from the supernova in the last spectrum of the 
data set ( day 1 248 ) . By comparing the WE data from day 1 to 
800 with the HST data (§ 7) from day 1278 to 2431, we esti- 
mated that at day 1248, only about 40% of our “background” 
subtracted SWP flux and about 57% of the LWP flux are in- 
deed from the supernova. Thus we estimated that 46% of the 
total “background”-subtracted WE flux at day 1 248 is from 
the supernova itself by adding the SWP and LWP flux compo- 
nents. Using this estimate and the observed total WE measure- 
ments from day 1 to day 1000 (excluding the noisy data points 
around day 940), we can fit a spline and estimate the UV con- 
tributions to the total UVOIR bolometric luminosity at epochs 
with OIR measurements. Notice that extrapolation of UV light 
curve beyond the last day of observation is required for the last 
CTIO and ESO data points (day 1 352 and 1316, respectively) . 
However, since the UV contribution to the total bolometric 
flux is less than 10% at late times, we estimate that the errors 
quoted in the Suntzeff et al. ( 1991 ) and Bouchet et al. ( 1991 ) 
papers will not be altered by the inclusion of the extrapolated 
UV fluxes. The extrapolated UV flux values for day 1316 and 
1352 are marked by crosses in Figure 17. The fitted spline to 
the UV light curve with the adjusted day 1 248 flux is shown in 
Figure 17 and the interpolated UV fluxes at dates with OIR 
measurements are presented in the second column of Table 4. 
The adjusted UVOIR bolometric luminosity data are pre- 
sented in the third column of Table 4A and 4B, respectively, 
for CTIO and ESO measurements. The inclusion of UV flux 
added less than 0.001 dex to the total bolometric luminosity 
by day 164. However, the contribution of UV increases with 
time — by day 468, inclusion of UV flux added 0.01 dex to the 
total bolometric flux. At the final date of calculation, day 1 352, 
the UV flux added 0.04 dex to the previous bolometric lumi- 
nosity calculations. 

As seen in Figure 1 7, both the bolometric luminosity and the 
U V flux from SN 1 987 A show an early decrease in magnitude. 
While the bolometric light curve reached a brief minimum at 
—day 8, the UV fight curve did not reach its minimum until 
—day 18 and it remained at the minimum for another month 
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before it started to increase. Despite differences in their prior 
developments, both light curves reached maxima at —day 90 
before they decreased again. By —day 1 20, the rate of decrease 
for UVOIR flux changed to the one suggested by the radioac- 
tive decay of 56 Co. The UV flux instead increased slowly to a 
maximum at —day 350 before it started to drop exponentially 
around day 600. Beyond —day 900, both curves seemed to 
have s slowed their rate of decrease. From Figure 17, it is obvi- 
ous that the UV fraction of the total UVOIR luminosity in- 
creases with time. The time development of this fraction is pre- 
sented in Figure 19. At the first day of WE observation (day 
1.66), emission in UV represents 46% of the radiation from 
SN 1 987 A. This fraction decreased steeply in the next few days. 
By day 1 8, less than 0. 1 % of the bolometric radiation was from 
the WE bands. The /u V /^uvoir ratio reached a minimal value 
of0.004% at —day 45. After that date, the contribution of UV 
flux started and continued to increase afterward. By day 800, 
UV accounted for —7% of the total UVOIR bolometric lumi- 
nosity of SN 1987A, and the compilation here is the only 
source where it has been included. The correct determination 
of the late-time bolometric flux from the supernova is espe- 
cially important to the identification its energy source. Frans- 
son & Kozma ( 1993 ) suggest that after —day 800, time-depen- 
dent effects due to long recombination and cooling times 
would lead the supernova ejecta to a “frozen” state, causing 
the emitted flux to be larger than the instantaneous radioactive 
energy deposition. Thus the observed supernova luminosity 
will be higher than that of the instantaneous energy input. 
Quantitative tests of this theory demand a set of correctly de- 
termined bolometric luminosity data, as given in Table 4. 

6. CONNECTION WITH DUST CONDENSATION 
6. 1 . Dust Condensation in SN 1 987 A Ejecta 

By comparing the atmospheric conditions in SN 1 987 A with 
classical novae, Gehrz & Ney ( 1987 ) suggested that dust grains 
would condense in the expanding ejecta when the temperature 
cools down to 1000 K. They predicted that the dust condensa- 
tion would occur by late 1987 or early 1988 (—day 300). A 



Days after outburst 

Fig. 19. — Fraction of UV emission to the total emission for SN 1987A. 
The integrated UV flux is compared with the adjusted UVOIR bolometric 
luminosity (Table 4) measured by CTIO ( filled circles) and ESO {open 
circles). 


similar study by Dwek ( 1988) made the same prediction and 
suggested that dust would form in the ejecta by t « 400 days. 
Hints of dust condensation came from the optical/ IR photo- 
metric and spectroscopic observations of the supernova. In 
Roche et al. ( 1989), they reported that the mid-IR continuum 
emission (8-13 pm), which had been fading steadily since 
—day 120, began to increase on day 465 and continued to do 
so up to day 578. However, the authors suggested that this 
emission originated from an IR echo. With subsequent obser- 
vations Roche, Aitken, & Smith ( 1993) concluded that the 
bulk of the mid-IR continuum flux originated instead from the 
newly formed dust in the ejecta which absorbed the optical 
photons and reradiated that energy in the IR region. Meikle et 
al. ( 1993) studied the 4 pm flux continuum and also suggested 
that dust condensation occurred as early as —day 350. At 
about the same time, the optical photometry data were fading 
at a rate faster (see Fig. 16) than the previous epoch. The in- 
crease in IR luminosity compensated for the optical deficit so 
that the overall bolometric light curve continued to follow the 
decay of 56 Co (Whitelock et al. 1989; Suntzeff & Bouchet 
1990). 

However, it was the spectroscopic evidence that was decisive 
in confirming the presence of condensed dust in the SN 1987 A 
ejecta. Around day 530, the line profiles of optical emission 
lines, such as Mg l] 4571, [O l] 6300, 6363, Ha, and [C I] 
9823, 9849, became asymmetrical and blueshifted by 500-700 
km s _I (Danziger et al. 1989). Lucy et al. ( 1989) explained 
this shift as the consequence of obstruction of the far side of 
the supernova ejecta by dust formed in the inner core of the 
supernova atmosphere, within the line-forming region. More- 
over, the bluer lines showed larger blueshifts, suggesting that 
some dust is optically thin. However, the condensation effi- 
ciency would be very small (^10‘ 3 ) if the only dust formed 
were small-sized grains (—0.01 pm) uniformly distributed in 
the ejecta. Lucy et al. ( 1991 ) resolved the problem by con- 
structing a model with dust grains concentrated in optically 
thick clumps with filling factor ^ 1 , while a diffuse distribution 
of small dust grains exists between these opaque clouds. They 
are able to obtain a better fit to the IR spectrum than their 
previous model (Lucy et al. 1989) in which only uniformly 
distributed dust grains are condensed. This clumping of the 
condensed dust is probably the reason the predicted optical- 
UV photometric “blackout” in Gehrz & Ney (1987) and 
Dwek ( 1 988 ) did not occur. This picture is consistent with the 
model constructed by Wooden et al. ( 1993) to match the IR 
observations from the Kuiper Airborne Observatory (KAO) 
from day 60 to day 775. In their model, the dust grains are 
clumped and these dust clumps emit a single-temperature 
gray-body spectrum in the infrared continuum of SN 1987 A 
so that most of the energy emitted redward of 6 comes from 

the dust. The dust temperature dropped from T g » 400 K at 
day 6 1 5 to T g « 265 K at day 775. Bouchet et al. ( 1 99 1 ) also 
fit a blackbody to the infrared photometric data for day 1316 
and obtained T g » 140 K. 

6.2. UV Color Evolution 

As seen in the optical + UV spectra in Figure 8, and the ratio 
of UV contribution to total bolometric luminosity in Figure 
19, SN 1987A is getting bluer in time. To investigate this 
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change quantitatively and to evaluate the effect of the forming 
dust on the UV emission, we define an ultraviolet “color” by 
subtracting the logarithm of integrated flux in the wavelength 
band 1600-2500 A from the logarithm of integrated flux in the 
band 2500-3300 A. The spectra in the range 1 1 50- 1 600 A are 
not included in the study because of their low signal to noise 
and the contamination by background stars. The time evolu- 
tion of this U V color is shown in Figure 20. After climbing to a 
early maximum at day 6 and stayed there until day 75, the UV 
color dropped, approximately in linear fashion, to a value of 
~0.8 by day 125. Then the color remained at about the same 
level for about a year (day 125-460) before it took a sudden 
turn toward the blue and started to decrease, more or less lin- 
early, until the end of our observing period. To examine the 
origin of the latter dip in color, the spectra just before and after 
day 460 are compared. The average spectrum for the period 
day 500-550 is divided by the average spectrum of day 350- 
400 and the result is presented in Figure 2 1 . The fluctuation of 
the ratio seen in the wavelength range 2200-2390 A is caused 
by the very low signal observed. Unfortunately, this is where 
part of the 2175 A “bump” of interstellar extinction is located. 
This bump, though not as prominent in the LMC as in the 
Galaxy (Fitzpatrick 1985), is the most prominent feature of 
the interstellar reddening function in the U V. However, by no- 
ticing the relative flux change in the 2040-2200 A range where 
stronger signals from the supernova are detected, we conclude 
that the 2175 A bump is not the dominant factor in the varia- 
tion of the U V color. In fact, the change of color is not caused 



Fig. 20. — The evolution of the UV color with log !0 ( time) (see text for 
definition of the UV color). Note that when the supernova gets redder in 
UV, it will have a higher value of UV color. 



Wav*4*ng1ti (Angstroms) 


Fig. 21. — Ratio of the average UV spectrum taken just before dust for- 
mation (day 350-400) over the average U V spectrum taken just after dust 
formation (day 500-550). (We assume that dust formation took place 
around day 460. ) The fluctuations in the range 2200-2 390 A are caused by 
the very low signal from the supernova in these two epochs. 


by the change of flux of any particular feature. Instead, it is 
caused by an overall change that seems to increase the flux to- 
ward shorter wavelength, resulting in a general “bluing” of the 
spectrum as time goes by. 

The time at which the U V color took a sharp turn to the blue 
is also the time when dust formation is believed to occur in the 
ejecta (day 450). Ordinary interstellar dust causes flux extinc- 
tion that increases with decreasing wavelength in the ultravio- 
let range, except for the distinct 2175 A bump. The fact that 
the supernova was getting bluer despite the dust formation sug- 
gests that the short-wavelength ultraviolet radiation did not 
suffer any extra extinction by the newly formed dust. In any 
case, the dust does not resemble interstellar dust in its extinc- 
tion properties in the UV. This is similar to what was seen in 
Nova Cassiopeiae 1993 in the UV where the formation of 
“gray” dust was observed (Shore et al. 1994). The solution of 
this problem lies in the composition, size, and distribution of 
the condensed dust. 

6.3. On the Reddening Properties of Dust in SN 1 987 A 

The supernova ejecta consist of regions of different elemen- 
tal compositions and physical properties, which might lead to 
the condensation of different types of grains. Lucy et al. ( 1 99 1 ) 
suggested that the condensation of “astronomical silicate” 
( Draine & Lee 1984) of radius 100 A, mostly in clumps, could 
explain both the steepening in the rate of decline of the [Si i] 

1 .65 ^m line flux and provide a better fit to the 10-20 fim con- 
tinuum. However, the [Si i] 1.65 pm line is very sensitive to 
changes in temperature (Fransson & Chevalier 1987), and 
thus its rapid decline may not be totally attributed to the con- 
densation to silicates. Using the theory of homogeneous nucle- 
ation and grain growth, Kozasa, Hasegawa, & Nomoto ( 1989, 
1991) modeled the condensation of clumped silicate and oxide 
grains (radius — 10-100 A), such as MgSi0 3 , A1 2 0 3 , and 
Fe 3 0 4 , to obtain fits to the observed 10 pm IR light curve. In 
the model by Wooden et al. ( 1993), the graybody radiation of 
the supernova in the mid- and far-infrared after ~day 400 is 
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explained by the presence of some heavily clumped dust clouds 
which are optically thick, with a total surface area filling factor 
of ~0.5. Spyromilio & Graham ( 1 992 ) suggested that the fall- 
ing of the gaseous Fe li near-IR emission could be explained by 
the combined effects of the rapid drop of temperature (below 
500 K), and that a fraction of Fe + got tied up in optically thick 
clumps of dust to escape the rapid cooling. In their efforts to 
model the light curves of neutral and singly ionized Fe, Ni, and 
Co infrared emission lines, Li, McCray, & Sunyaev (1993) also 
constructed a model of the supernova atmosphere with some 
iron-rich clumps with a total filling factor of ~0.5. This sug- 
gests that the newly formed dust may be rich in iron-group 
elements. This conclusion is supported by the study of Wick- 
ramasinghe & Wickramasinghe ( 1993). They argued that in 
the heavily non-LTE conditions of the supernova atmosphere, 
the condensation of metallic iron might precede the formation 
of iron oxide, thus leading to the production of “iron whis- 
kers.” These metallic iron needles of radius 100 A and length 1 
mm are shown to have the necessary mass absorption coeffi- 
cients to explain the lack of observed flux in the far-infrared to 
microwave range ( 3 /xm to 30 cm) in the day 1 300 spectrum of 
SN 1987A. As described in § 3, the UV spectrum ofSN 1987A 
is dominated by the heavily blanketed lines of iron group ele- 
ments, especially Fe n. The combined effects of decreased 
opacities of iron lines with the formation of iron rich dust, and 
increased absorption by the condensed dust, may explain the 
evolution of the UV color near the epoch of dust formation. 

As described in § 3, the UV radiation from the supernova 
went through a series of resonance scatterings before escaping 
from the envelope. Fransson ( 1 994 ) pointed out that since the 
effects of dust absorption are prominent for optical lines, we 
should see a larger fraction of the UV radiation, with a much 
larger path length due to resonance scattering, to be absorbed, 
thermalized, and emerge in the infrared after dust condensa- 
tion. However, such a large drop in flux in UV is not observed 
(see Fig. 16). In fact, the supernova continued to grow bluer in 
the U V even as dust condensation proceeded. This effect might 
be caused by dust grains with decreasing absorption coeffi- 
cients with decreasing wavelength in the U V. For small grains 
(radius £ X / 10), Rayleigh scattering is the dominant scatter- 
ing process and its effect goes with the inverse fourth power 
of the wavelength. We found that neither the “astronomical 
silicate” ( radius — 1 00 A ) proposed by Lucy et al. (1991), nor 
the oxides and silicates grains (radius ~~ 10-100 A) proposed 
by Kozasa et al. ( 1991 ), nor the iron whiskers (radius 100 A; 
length 1 mm) proposed by Wickramasinghe & Wickramasin- 
ghe ( 1993), has the correct absorption cross section depen- 
dence with wavelength to explain the bluing effects in the UV. 
Thus we would need different components in the dust other 
than the these small sized grains. With larger sized grains, Mie 
scattering will be increasingly more important. If the dust is in 
front of the UV source, the effect of forward scattering could 
work to make the spectrum bluer. Lucy et al. ( 1 989 ) derived a 
constraint on the size of the dust grains (“astronomical sili- 
cate” radius ^ 500 A) to explain the bluing of the optical emis- 
sion lines. We hope a constraint on the grain size can be deter- 
mined with the UV color information. 

One alternate explanation for the bluing of the spectra in 
UV following dust formation is as follows: dust condensation 
started in the cooler part of the ejecta. If some of the dust 


clumps formed are totally black, the newly condensed dust will 
effectively block out the lowest temperature material in the su- 
pernova envelope. Therefore, even though the flux from the 
supernova was declining, the temperature of the UV radiation 
could go up. And as a consequence, the UV color would get 
bluer in time. A detailed inhomogeneous model of the super- 
nova debris is needed to see whether this effect can account for 
the observations. 

1. COMPARISON WITH HST RESULTS 

The IVE data form an essential record that will prove useful 
for interpreting future observations of SN 1 98 7 A with the Hub- 
ble Space Telescope (HST). Observation of SN 1987 A was 
part of the science assessment program to test the performance 
of the telescope. Continuing observations of the supernova 
were carried out by the Supernova INtensi ve Study ( SINS ) GO 
program. The Faint Object Camera (FOC) aboard the HST 
has the resolution capability (0T022) to resolve the supernova 
from its neighboring stars. It also has the sensitivity to detect 
the fading supernova after it has become too faint for IUE mea- 
surements. The first HSTUV observations of SN 1987 A were 
taken on 1990 August 24 (day 1277.8) with the FOC f/96 
through the F175W and F275W filters ( Jakobsen et al. 1991 ). 
Further results on the expansion of the supernova debris and 
photometry were also published (Jakobsen, Macchetto, & Pa- 
nagia 1993) based on a second set of observations made at 
1991 December 13 (day 1754.3), and also in Jakobsen et al. 

( 1994) based on the first set of post-COSTAR observations. 
Complete photometric results from the SINS SN 1 987A obser- 
vations will be presented in Kirshner et al. ( 1995 ) . 

The UV photometry in the HST bands for epochs before 
HST deployment can be obtained by convolving the IUE spec- 
tra with the FOC filter passbands. We convolved the IUE spec- 
tra with the FOC f/ 96 F 1 7 5 W filter function ( peak wavelength 
1730 A, FWHM 714 A) and F275W filter function (peak 
wavelength 2740 A, FWHM 594 A) using the CALCPHOT 
task in the SYNPHOT program of the Space Telescope Science 
Data Analysis System (STSDAS) (Bazell 1990). The CALC- 
PHOT task is capable of synthesizing FOC photometry corre- 
sponding to the selected HST filter bands with the input spec- 
trum using equation ( 1 ). The input UV spectra are corrected 
for neither the parent galaxy redshift nor interstellar extinction 
so that the derived spectrophotometric data can be compared 
directly to the photometric measurements. The synthesized 
FOC photometry results for the filters F175W and F275W are 
shown in Figures 22a and 22 b, respectively. 

The synthesized photometry results from IUE spectra can 
be compared with the photometry directly measured from the 
HST FOC images obtained under the SINS program. It is for- 
tunate that the first UV observation of SN 1987A by HST was 
made not long after the signal from the supernova in UV was 
too weak for detection by IUE. Six sets of observations have 
been made at day 1278, 1754, 1876,2031,2260, and 2431 with 
images from both the F175W and F275W filters (except at day 
1754, when only the F275 W image is taken ) . The complete log 
of the FOC UV data set is given in Table 5. In all the FOC 
images, the SN 1 987A is superimposed on the outer part of the 
point spread function (psf) of star 2 and star 3, which are both 
well saturated in the center in all the frames. The flux from the 
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Fig. 22. {a)-(b) The FI 75 Wand F275W light curves of SN 1987A from the IUE and HST observations. The filled squares are the results from the 
SYNPHOT program by convolving the I UE data with the HST observing passbands. The open circles are the SYNPHOT results of convolving the combined 
UV + optical data with the HST passbands. The filled circles represent the measured photometry from HST images. 


two asymmetric halos represents the major source of contami- 
nation to the correct photometric measurement for the super- 
nova. The DAOPHOT II program was used to create a model 
of the entire V psf of star 2 (Stetson 1991 ). This model was 
subtracted from the HST image. The same model was scaled 
to the flux magnitude of star 3 and subtracted from the image. 
This procedure is valid despite the fact that both star 2 and 
star 3 are saturated at the center of the FOC images for both 
passbands because we are only interested in determining the 
flux from the supernova, which is located at the periphery of 
the psf of the two stars. The flux from the supernova is mea- 
sured by integrating over a small circle (5 pixels, OH l ). The 
result is multiplied by the fraction of encircled energy (22%) in 
the psf at the radius of 0T1 1 . This fraction is determined from 
the high signal to noise images of SN 1992 A (Kirshner et al. 
1993) obtained with the FOC f/96 F175W and F275W filters. 
These images provided good measurements of the encircled 
energy as a function of psf radius. The absolute UV photome- 
try measurements by FOC is correct to — 10%- 1 5% ( P. Hodge, 


TABLE 5 

HSTUV Observations of SN 1987A with 
Faint Object Camera (FOC) 


Date (UT) Day Filter (s) 


1990 Aug 24.0 1277.7 F175W 838 + 838 

1990 Aug 24.1 1277.8 F275W 838 

1991 Dec 13.6 1754.3 F275W 1200 + 838 

1992 Apr 13.7 1876.4 F175W 2000 

1992 Apr 13.7 1876.4 F275W 2000 

1992 Sep 15.7 2031.4 F175W 2000 

1992 Sep 15.7 2031.4 F275W 2000 

1993 May 1.9 2259.6 F175W 2500 

1993 May 1.8 2259.5 F275W 2500 

1993 Oct 20.3 2431.0 F175W 1500+ 1700 

1993 Oct 20.3 2431.0 F275W 1450+1550 


private communications) with the major uncertainty coming 
from the determination of the standard star flux. The photom- 
etry measurements in this present work are — 1 5%-40% lower 
than the previously published results from the day 1278 and 
day 1 7 54 data ( Jakobsen et al. 1991, 1 993 ) . The discrepancy 
is probably caused by the different methods employed in re- 
moving the fluxes from star 2 and star 3 in the images. 

There is one detail not to be overlooked while comparing 
the IUE and HST results— both the FI 75 W and F275W filters 
have extensive red leaks that include optical light beyond the 
observed limit of the IUE cameras at 3300 A (Nota, Jedrzej- 
ewski, & Hack 1993). In Figure 22, the open circles represent 
the synthesized FOC photometry results obtained with the 
SYNPHOT program with IUE data ( 1 150-3300 A) as the in- 
put. To study the effects of the red leaks, the results from the 
SYNPHOT program with the combined UV 4- optical data 
( 1 150-1 1000 A) as the input spectra are also plotted (filled 
squares). We can see that neglecting the supernova flux long- 
ward of 3300 A can cause a significant underestimate of the 
UV photometry. At day 44, excluding flux above 3300 A will 
result in FOC-band UV photometry 35 and 3 times lower than 
the full photometry for F175W and F275W filters, respec- 
tively. For the FI 75 W filter, even though the transmittance be- 
yond 4000 A is less than 0.3% of that at the transmission max- 
imum at 1730 A, on day 44 the flux in the optical regime is 

1 0 4 times larger than that at 1 7 30 A. The overall effect is that 
the UV photometry measurement with the FI 75W filter at day 
44 would have about 97% of the flux originated redward of 
3300 A. Similarly, for the F275W filter, while the transmit- 
tance beyond 4000 A is less than 0.75% of that at 2740 A, the 
flux in the optical regime is about 500 times larger than that 
around 2740 A at day 44. Thus only about one-third of the 
collected light through the F275W filter would have come from 
the UV. As mentioned in § 5.2, as the supernova ages, rela- 
tively more flux of SN 1 987A comes from the ultraviolet region 
as compared to optical. Therefore, the difference between the 
SYNPHOT results using UV 4 optical input and UV only in- 
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put diminishes with time. For observations after '—day 700 for 
F175W and after ^day 400 for F275W, the underestimate of 
synthesized photometry caused by the cutoff at 3300 A be- 
comes negligible. 

The WE and HST results agree well for measurements 
made with the F275W filter around day 1250, the time at 
which the WE satellite was not sensitive enough for further 
detection of the dimming SN 1987A, and when we started our 
program with HST on SN 1987 A. This confirms that our 
method of background stellar subtraction for WE is accurate, 
as compared to the one employed by SFdC93, which errone- 
ously generated a late-time flux — 10-15 times higher than our 
study. Our WE and HST F175W results agree only moder- 
ately well around day 1 250, with the WE measurement about 
0.5 dex higher than the HST results. This is probably caused 
by the additional flux detected from the variable star 3 as de- 
scribed in § 5.2. 

8. CONCLUSION 

We have described the details to create a complete set of ul- 
traviolet spectra taken by the WE satellite on SN 1987A. We 
hope that this data set will contribute toward understanding 
the late stage of stellar evolution, the supernova explosion 
mechanism, and evolution of the explosion debris. The major 
problem we faced in the data reduction process has been the 
removal of the flux contributions from nearby field stars. No 
previous analysis has had background measurements of suffi- 
cient accuracy. Our WE results are consistent with other UV 
observations such as Walraven VBLUW photometry and HST 
images. 

The complexity of the line blending in the U V region implies 
that detailed analysis is needed to obtain abundances and even 
line identification in the supernova atmosphere. Early WE 
and optical spectra had been used for modeling the SN 1987A 
atmosphere (Harkness & Wheeler 1988; Eastman & Kirshner 
1 989; Schmutz et al. 1 990; Hauschildt & Ensman 1 994 ) to un- 
derstand the effects of scattering and geometry of the expand- 
ing envelope. LTE synthetic spectral analysis for the optical 
and UV spectral data taken by HST had been done on SN 
1992A, a SN la (Kirshner et al. 1993), and SN 1993J, aSNII 
(Jeffery et al. 1 994 ) . The U V behavior of supemo vae is closely 
related to the degree of interaction of the supernova ejecta with 
the circumstellar materials. For supemovae with thick circum- 
stellar envelopes, such as SN 1979C, SN 1980K, and SN 
1993J, relatively smooth UV spectra are observed, whereas for 
SNe la and SN 1987A, which had a weaker radio display and 
hence a more dilute circumstellar medium, showed UV-defi- 
cient spectra with broad band absorption caused mostly by sin- 
gly ionized iron group elements. The UV deficit is caused by 


expanding atmospheres with moderately steep density profiles, 
together with a huge number of resonance lines of the iron 
group elements providing the U V opacity. These closely sepa- 
rated lines blended under the expansion and caused the line- 
blanketing effect, and the broad features observed will natu- 
rally emerge from the blending of the strongest absorption 
lines. 

With the WE data now swamped by star 2 and 3 (both B 
stars), future long-term monitoring of SN I987A in UV de- 
pends on observing programs with the HST The HST correc- 
tive optics will provide significant improvements in photome- 
try by confining the flux from star 2 and star 3 to a small region 
far from SN 1 987A. While post-COSTAR FOC or WFPC2 im- 
ages can provide continuous UV light curve coverage 
( Jakobsen et al. 1994), the FOS should provide spectroscopic 
data on both the expanding debris and the circumstellar ring 
(Wang et al. 1994). Even though the supernova is fading in 
UV at present, it can become bright again in the future. The 
opaque UV atmosphere may turn transparent and reveal the 
stellar interior composition of the progenitor, Sk -69°202, 
right before the explosion. The UV spectrum then becomes a 
emission-line spectrum, and it should prove to be helpful in 
determining the masses of newly produced elements such as C, 
Mg, and Si, which all have strong lines in the U V region. 

The circumstellar ring can also become a bright source of 
UV radiation when it is hit by the envelope of SN 1987 A. Luo, 
McCray, & Slavin ( 1 994 ) predicted that this will occur in 1 2 ± 
3 years after explosioirfAD 1999 ±3). The shocked ring will 
become a bright source of optical, UV, and X-ray emission 
(Luo Sc McCray 1991; Luo et al. 1994; Suzuki, Shigeyama, & 
Nomoto 1993; Masai & Nomoto 1994). Ultraviolet emission 
lines of He, C, N, and O will all be rising in brightness ( up to 
— 10" 12 ergs cm -2 s" 1 ) when the shock occurs and will stay 
bright a few year after the shock (Luo et al. 1994). The time 
evolution of these emission lines can be studied in detail with 
future HST observations. 

We are extremely grateful to the staff at both the GSFC and 
VILSPA WE. Observatories for supporting the continuous 
observing program of SN 1 987 A by the WE satellite, especially 
Yoji Kondo for his energetic advocacy of this work. We would 
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3; J. W. Pel for the VBLUW photometry; L. S. Lyubimkov for 
the ASTRON spectra. C. S. J. Pun would also like to thank D. 
Jeffery, B. Schmidt, and P. Hoflich for many insightful discus- 
sions on SN 1987A; T. Kozasa, and N. C. Wickramasinghe for 
their comments on dust condensation; S. M. Dougherty for the 
discussions on Be stars. This work was supported in part by 
NASA grant NAG 5-841. 


APPENDIX 

We have included the full set of processed WE spectra on SN 1987A (both LWP and SWP) listed in Tables 2 and 3 in the AAS 
CD-ROM Series, Volume 5. The combined optical T UV spectra presented in Figure 8 are also included. Please consult the 
documentation included in the data file for further information. 
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